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LIST OF SYMBOLS 



g 



n 



of 



A 



A 



w 



a 



Area of a structural member cross section (sq. in.) 
Tributary area of a structure for purpose of computing load 
(sq. ft) 

Area of steel beaan cross section, composite beam construction only 
(sq. in.) 

Net area of the rear face of a rectangular structure with openings 
(sq. ft) 

Area of concrete in cross section (sq. in.) 

Area of flange of structural steel beam or girder (sq. in.) 
Net area of the exterior front wall of a rectangular structure 
with openings (sq. ft) 

Cross section area of concrete column (sq. in.) 

Gross area of the exterior front face of a rectangular structure 
(sq. ft) 

Area of each portion of the subdivided face of a rectangular struc- 
ture with openings (sq. ft) 

Area of openings in the front wall of a rectangular structure 
(sq. ft) 

Net area of the back wall of a rectangular structure with openings 
(sq. ft) 

Area of reinforcing steel bars parallel to steel beams in composite 
construction (sq. in.) 

Area of tension steel in reinforced concrete member (sq. in.) 

Area of one stirrup of the web reinforcement of a reinforced .con- 
crete member (sq. in.) 

Area of web of steel beam (sq. in.) 

Area of compression steel in reinforced concrete member (sq. in.) 

Clear distance between flanges of a steel beam (in.) 
Depth of compression area of concrete fin.) 
Short side dimension of two-way slabs (ft) 
Column spacing in flat slab construction (ft) 
Acceleration of mass (ft/sec^) 

Acceleration of the g floor mass at times t.,. 



respectively 



th 



a (t) Acceleration of the g floor mass as a function of time 



a ~ . ... a 



.., 
n- 



n 



rn-l 



Acceleration at time 
respec tively 
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a 1 Distance from edge !t a !t to centroid of the load on area "A" m 
a" Distance from edge "a" to the centroid of the inertial force 

B Peak value of externally applied load 

Elastic stiffness coefficient for soil under rocking foundation 

B Equivalent peak value of the externally applied load 




B Peak value of the effective load on the g th floor mass of a 



multi-story building 
B Peak value of the external load, f (t), on the g floor mass 

o g 

b Width of flange of steel beam (in.) 

Width of compression flange of reinforced concrete beam (in.) 
Long side dimension of two-way slabs (ft) 

b f Flange width of channel section (in.) 

b Width of longitudinal web stiff ener for steel beams (in.) 

b f Width of stem of reinforced concrete tee beam (in.) 

Width of flange of steel beam used in composite construction (in.) 
Distance from edge "b" to centroid of load on area "B" 

b" Distance from edge "b" to the centroid of the inertial force 

C Compression force developed in reinforced concrete member (kips) 
Static compression column influence factor for shear walls 

C Crater depth factor for underground burst 

C, Dynamic compression column influence factor for shear walls 
Average drag coefficient 

C Local drag coefficient 

C Maximum ACI Code moment coefficient obtainable for any point on 
a slab (equation 6.52) 

C R Ratio of maximum resistance to peak load, C = R /B 
C R Parameter defined by equation 59^ 

C Seismic wave velocity (fps) 
s 

C T Ratio of load duration to natural period of oscillation, C T = T/T 

C, Approximate coefficient for determining elastic truss deflection 

(varies with type of truss) 

w 

C TT Ratio of maximum work done to absolute maximum work done, C TT = ri 
W W W 

P 
C' Local drag coefficient for a cylindrical segment of infinite length 

C. . Value of determinate formed by omitting the row and column con- 

1J taining the coefficient C. . from the matrix 

^- J 
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c Distance from neutral axis to the extreme fiber (in.) 

Cohesive strength of the soil as determined by conventional 

laboratory testing procedures (psi) 

The clear distance between compression flange and stiffener 

c Velocity of sound in undisturbed air (fps) 

c f 1 Velocity of sound in the region of reflected overpressure (fps) 

D Diameter of a cylindrical segment or a spherical dome (ft) 
Diameter of the smallest longitudinal reinforcing bar (in.) 
Diameter of spiral core of a concrete column 

D.L.F Dynamic Load Factor = x /x 

m s 

D.I.F Dynamic Increase Factor for strength of materials 

d Depth of burst of bomb below ground surface (ft) 
Height of burst of bomb above ground surface (ft) 
Depth from the compression face of a beam or slab to the centroid 
of longitudinal tensile reinforcement (in.) 
Least lateral dimension of rectangular column (in.) 
Total depth of steel beam (in.) 
Over -all depth of truss (ft) 

Diameter or width of column capital in flat slab construction (ft) 
Diameter of circular reinforced concrete column (in.) 

d 1 . Distance between centroids of compression and tensile steel in 
doubly reinforced concrete member 

E Modulus of elasticity, including Young's Modulus for steel and 
soil, of any structural material except concrete 
Energy absorbed by the equivalent system 
Explosive factor for cratering 
Compressive modulus of elasticity of concrete, kips/sq. in. 

E Initial tangent modulus of elasticity for concrete 
E Explosive factor for underground burst 

"t~V 

E Maximum allowable energy absorption for the g story 

o 

e Base of natural logarithms (2.718) 

Eccentricity of load, distance from gravity axis to point of 

application of load (in.) 

Center to center spacing of two trusses 

e Strain at failure (%) 

c , 

e Strain at initiation of strain hardening ($) 
s 

e Strain at yield ($) 
j 

F Horizontal force 

F /,\ Sum of all external loads above the g story, f(t) + f (t) +..., 
8( ' + f (t) S g 

F Design lateral load on frame at time t 

IX 
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F Q Summation of all external horizontal forces applied to structure, 
including foundation reactions 

F Pressure factor for underground burst 

Normal component of total passive resistance force (kips) 

F(t) Load applied to a structural element or system as function of time 

f Fiber stress (psi) 

f a Column critical buckling stress (psi) 

f ad Average stress at design level due to axial load (psi) 

f B Dynamic yield strength for structural steel in composite beam 
construction only (psi) 

f^ Elastic buckling stress due to bending (psi) 

f t(i Maximum stress due to bending at design level for steel members 
subjected to axial load and bending (psi) 

P, Dynamic yield strength of steel (psi) 

r gx< 



f Y (t) Effective load on the g floor of a multi-story structure as a 



function of time 
f.^ Initial static stress (psi) 
f lv Lover static stress (psi) 

f Dynamic yield strength for reinforcing steel for composite beam 
construction only (psi) 

f , Tensile strength of concrete (psi) 

f(t) Load applied to a structural element of system as a function of 
time 

f Upper static yield point stress (psi) 

f Static yield point stress 
j 

f _ , f p , f Factors used in the design of tee beams 

3 4-V> 

f-^t), f 2 (t), f (t) Load applied to the first, second, and g floors 

of a multi-story building as a function of time 

f f Static ultimate compressive strength of concrete 
f ' Dynamic ultimate corapressive strength of concrete 

G Solidity ratio of truss 

Shear modulus of elasticity (psi) 

G.Z. Ground zero (point on ground below point of detonation of air 
burst) 

g Acceleration of gravity 

H Impulse per unit area, impulse per foot of length, or total impulse 
(kip-sec) 

Depth of footing or vail below the surface of the r^ound 
Height of shear wall to center of roof beam 
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H Equivalent impulse acting on the equivalent system (kip-sec) 

f-Vt 

H . Total impulse of the external load on the g floor mass (kip-sec) 
H, Total impulse of the external load (kip-sec) 

u 

h Height (ft) 

Story height (ft) 

Column height (ft) 

Height of a structure above ground (ft) 

Maximum thickness of channel flange (in.) 

h Clear height of column (ft) 

*f"V* 

h Height of the g story of a multi-story building (ft) 

o 

h Clearing dimension for each portion of the subdivided front face 
n of rectangular structure with openings (ft) 

h 1 Clearing height of a closed rectangular structure (ft) 

Vertical distance from roof to point under consideration for buried 
structure (ft) 

hJ Weighted average buildup height of the exterior rear face of a 
rectangular structure with openings (ft) 

h' Weighted average clearing height for the exterior front face of 
a rectangular structure with openings (ft) 

h' Weighted average clearing height for the interior rear face of a 
rectangular structure with the openings (ft) 

h! f Weighted average buildup height for the interior front face of a 
rectangular structure with openings ('ft) 



I k h\ 
Moment of inertia 1 in. or ft J 



I. Inertial force of "A" portion of two-way slab 

I Average of the gross and transformed moments of inertia 



(in. or ft ) 



3L Inertial force of "B" portion of two-way slab 

B 



I Moment of inertia of the gross section I in. or ft 

g f k k\ 

I. Moment of inertia of the transformed section ^in. or ft j 

I. Moment of inertia of rectangular section in concrete tee beam 

1 I M 

design I in. J 

IP Moment of inertia of tee section in concrete tee beam design 
I =* I ft ) Velocity impulse at time, t 

I Mass moment of inertia of structure about axis of rotation "0" 
o 

j Ratio of distance between centroid of compression and centre id 
of tension to the depth, d 

K Modulus of strain hardening 
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KE Kinetic energy ' 

(KE) Kinetic energy of the actual system 

a 

(KE) Kinetic energy of the equivalent system 

K Either load, mass, or load-mass factor for slab fixed on four sides 

K. ., (i = 1, 2, 3,...; j = 1, 2, 3 ) Stiffness influence coefficient 

1J (klps/ft) 

1C Load factor 

J/L Load-mass factor 

K, Mass factor 

Kp Normal component of passive pressure coefficient accounting for 

the cohesive effect of any soil having an internal friction equal 
to 

Kp Normal component of passive pressure coefficient for any soil with 
internal friction angle equal to and with zero wall friction 
developed 

Normal component of passive pressure coefficient for any soil with 
internal friction angle equal to and with maximum wall friction 
developed 

K- Resistance factor 

K Either load, mass, or load-mass factor for slab simply supported 
on four sides 

KT Kiloton, 1000 tons 

K Either load, mass, or load-mass factor for special edge conditions 
z 

K Q Ratio of the maximum average overpressure to the reflected over- 
pressure existing on an inclined roof 

K 1 Beam equivalent length coefficient used in the design of steel 
"beam colujftns 

K !l Column length factor used in the design of steel columns and "beam 
columns 

k Spring constant, force required to cause unit deflection of spring 
(kips/ft) 
Soil pressure factor (psi) 

k Soil constant for underground explosion (psi) 
k- Effective spring constant (kips/ft) 
k Equivalent spring constant (kips/ft) 

k Spring constant in the elasto-plastic range 

e P ^jj 

k-, k~, ...k ... Spring constants for the first, second, and the g 
g stories of a multi-story building (kips/ft) 

k Simplified notation of k 
6 
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Symbols 



. 

sl 

V 1 

gi 

k.., 




L f 



M 



M 
M 






M 



cp 



M- 



, 



, 



M 



m 



f 



M 

^ 



M 



mp 



-Spring constant of the coupling spring between the g 
floors of a multi-story "building (kips/ft) 
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and the i 



gi 



- V 



k , .k 



. 
g 



gi 
k 



Spring constants of the springs connecting m with 



the masses nu, nu, . . .m.,...m 






Span length of beam or truss (ft) 

Unsupported length of beam (ft) 

Length of* a rectangular structure in the direction of propagation 

of the blast wave (ft) 

Length of a cylindrical segment along the cylinder axis (ft) 

Length of shear wall center to center of column steel 

Spacing of columns in each direction 

Distance from the outside of the front face to the inside rear 
face of a rectangular structure with openings (ft) 

Distance from the front face of a rectangular structure to a point 
under consideration on the roof or sides in the direction of pro- 
pagation of the blast wave (ft) 

Length between sections of zero and maximum moment being considered 
(composite beams) 

Bending moment applied to a section 

Moment of forces on piles about their centroidal axis 

Resisting moment of soil on the footing per unit width 

Bending moment at center line of a beam or slab 

Negative resisting moment in column per foot, flat slabs (kip-ft) 

Positive resisting moment in column strip per foot, flat slabs 






Maximum design moment in a member under axial load, P 

Dynamic bending moment 

Elastic dynamic buckling moment 

Maximum bending moment 

Component in a plane perpendicular to edge "a" of the total re- 
sisting moment along the fracture lines bounding area "A," two- 
way slabs (kip-ft) 

Component in a plane perpendicular to edge "b, n of the total re- 
sisting moment along the fracture lines bounding area B, two-way 
slabs (kip-ft) 

Negative resisting moment in middle strip per foot, flat slabs 



Positive resisting moment in middle strip per foot, flat slabs 
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M Negative bending moment at support 

M Moment of all external forces on structure about axis of rotation 

O 

M f Summation of external moments about the point of rotation ex- 
cluding footing projection moments 

Mp Plastic resisting moment under bending only 

ML Negative plastic resisting moment in column strip per foot, flat 
Fcn slabs (kip-ft) 

Mp Positive plastic resisting moment in column strip per foot, flat 
p slabs (kip-ft) 

Component of the total plastic bending moment capacity along the 
fracture line boundary of area "A" which is in a plane perpendic- 
ular to edge "a," or the total positive plastic bending moment 
capacity for a section parallel to edge "a," two-way slabs (kip-ft) 

Component of the total plastic bending moment capacity along the 
fracture line boundary of area "B IT which is in a plane perpendic- 
ular to edge "b," or the total positive plastic bending moment 
capacity for a section parallel to edge "b," two-way slabs (kip/ft) 

HPm Plastic resisting moment at centerline of beam or slab 

M-, Negative plastic resisting moment in middle strip per foot, flat 
slabs (kip-ft) 

ML. Positive plastic resisting moment in middle strip per foot, flat 
^ slabs (kip-ft) 

M Maximum positive bending moment (kip-ft) 
pos 

Mp Plastic resisting moment at support (kip-ft) 

Mp Total negative plastic bending moment capacity along edge "a," 
a two-way slabs (kip-ft) 

Total negative plastic bending moment capacity along edge tf b, ff 
two-way slabs (kip-ft) 

Yield resisting moment (kip-ft) 

Plastic resisting moment per unit of width of slab (kip-ft/ft) 

Plastic positive bending moment capacity per unit width for short 
span, two-way slabs (kip-ft/ft) 

Plastic negative bending moment capacity per unit width at center 
sa of edge "a" for long span, two-way slabs (kip-ft/ft) 

Mp , Plastic negative bending moment capacity per unit width at center 
s of edge if b" for short span, two-way slabs (kip-ft/ft) 
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ML Moment at the intersection of the two-linear portions of the P vs M 
curve of a steel beam section (kip-ft) 

MI Theoretical plastic resisting moment (kip-ft) 

f P 2\ 

m Miass per unit length ^kip-sec /ft j 

Point mass ^kip-sec /ft] 

In concrete design: f /0.85 f ' 

Total moving mass of structure and earth included between footings 
Number of fundamental dimensional quantities 

/ 2 2\ 

m Mass of the equivalent system ^kip-sec /ft j 

m Mass of structure considered to rotate as well as translate 

m t Total mass of the element or structural system under consideration 

(kip-sec 2 /ft) 
m. Total mass of the equivalent system (kip-sec /ft) 

m., nu-.m Mass of the first, second, and g floors of a multi-story 
- 1 ^ g building 

N Number of non-dimensional parameters 

Total number of shear connectors required between the points of 
zero and maximum moment of a composite beam 
Weighted number of piles in group 

N f Ratio of the length of fixed-edge perimeter to the total perimeter 

N g Ratio of the length of simply- supported perimeter to the total 
perimeter 

n Number of columns in a story 

Ratio of modulus of elasticity of steel to modulus of elasticity 

of concrete 

Number of dimensional variables 

Number of stories in a multi-story frame 

P Load or force (kips) 

Total dynamic load on slab (kips) 
Panel influence factor for shear walls 

P A Total load on "A" portion of two-way slab (kips) 

P Average axial load acting on each of several columns of a frame 
(kips) 

P, k Local overpressure on the back face of a buried rectangular 
structure (psi) 

P B Total load on "B" portion of two-way slab (kips) 
P Compression mode overpressure (psi) 

P Uniform compressive pressure applied radially on arch 
ca 

P Uniform radial pressure that produces buckling of arch 

C A 
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P C Local overpressure normal to the exterior surface of a cylindrical 
segment (psi) 

P D Maximum axial load on column with given M^ (kips) 
P' Reduced maximum axial load for long columns 

P, Deflection mode overpressure (psi) 
Total dynamic axial load 

P dome Local overpressure normal to the exterior face of a spherical 
dome (psi) 

P e Equivalent concentrated load for equivalent system (kips) 

P , Elastic dynamic buckling load (kips) 

P (t) Equivalent load on an element as a function of time 

P f Average value of P(t) for the far (or leeward) side of the arch 
(psi) 

Local overpressure on the front face of a buried rectangular 

i / . \ 
structure (psi) 

P Peak underground overpressure resulting at a given location from 
^ an underground burst (psi) 

P Total vertical load (blast plus static) on column 
Force acting at any time, t 

Average value of P(t) for the near (or windward) side of the arch 
(psi) 

P Dynamic plastic axial load capacity of column 

P f>1 Reflected shock wave overpressure for angle of incidence of zero 
degrees (psi) 

P Reflected shock wave overpressure for angle of incidence other 
than zero degrees (psi) 

P Overpressure existing in the incident shock wave for any value of 
t-t d (psi) 

P , Overpressure existing in the incident shock wave when t-t , = t, 
Sb (psi) d ^ 

P . , Local overpressure on the exterior side walls of a rectangular 
structure 

P Initial peak incident overpressure (psi) 
so 

P . Peak overpressure of the shock wave formed in the interior of a 
structure with openings (psi) 

P , Stagnation overpressure; the overpressure existing in a region in 
^ which the moving air has been brought completely to rest (psi) 

P Axial load on the columns of a frame due to the vertical live and 
dead loads (kips) 
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The time average of the total vertical load on the column in the 

time interval t < t < t 

ge gm 

P Average overpressure on the exterior back wall of a rectangular 



, 
ack 



structure (psi) 



P, , Average net overpressure acting on back wall of a rectangular 
b ~ net structure (psi) 



( P back) 



Peak value of the average overpressure on the exterior back 
ma:x: wall of a rectangular structure (psi) 

Average overpressure on the closed ends of a cylindrical segment 



_ 

P . Average net overpressure acting on front wall of a rectangular 



' f -net 



structure (psi) 



P , Average overpressure on the exterior of the front wall of a 
front , - , , / . \ 
rectangular structure (psi) 

P- f + Average overpressure on the interior front wall of a rectangular 
structure with openings (psi) 

P, - Reflected shock wave average overpressure in the interior of a 
rectangular structure with openings (psi) 

P. f Average overpressure on the interior roof of a rectangular struc- 
ture with openings (psi) 

P. . , Average overpressure on the interior side walls of a rectangular 
structure with openings (psi) 

P . Net average horizontal overpressure exerted on a rectangular 
structure (psi) 

P Average net overpressure acting on the roof of a rectangular 
r net , , / . \ 
structure (psi; 

P Average overpressure on the front slope of a gable roof (psi) 

Average overpressure on the exterior of a rectangular structure 
(psi) 

P ? Peak average overpressure on the front slope of the first gable 
of a multi-gabled roof (psi) 

? , Average net overpressure acting inward on sidewall of rectangular 
structure (psi) 

P . , Average overpressure on the exterior side walls of a rectangular 
side . r / . \ 
structure (psi; 

P, Axial load determined by the intersection of the two linear 
portions of the P vs M curve of a steel beam section (kips) 

P(t) Actual load on a structural element as a function of time 

p Ratio of tensile reinforcement in reinforced concrete members to 

concrete area, A /bd 

Uniformly-distributed load intensity 
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p Critical steel ratio for reinforced concrete member 
p Equivalent static load for an arch (psi) 

p ! Ratio of volume of spiral reinforcement to the volume of the 

concrete core (out-to-out of spirals) of a spirally reinforced 

concrete column 

Ratio of compressive reinforcement in beams to the concrete area, 

A!/bd 
s 

Q Statical moment of the section about the centre idal axis 
Strength of a shear connector in composite construction 

/ o\ 

Radiant energy on a unit area I cal/cm j 
q. Unit drag pressure produced by the incident shock vave (psi) 

q. Dynamic design bearing pressure on soil ( kips/ft *~ ) 

/ V 2\ 
q Maximum bearing pressure on soil (kips/ft y 

q. Maximum unit drag pressure produced by the incident shock wave 
(psi) 

R Radius of a spherical dome 

Total resistance of structural element or structural system (kips) 
Dosage of gamma radiation (roentgens) 
Crater radius (ft) 

R. Total resistance of M A T! portion of two-way slabs (kips) 

R Crater radius (ft) 
u 

R Horizontal static load resistance of shear wall at first cracking 
R, Horizontal dynamic load resistance of shear wall at first cracking 

R Reynolds number of the high velocity wind in the incident shock 
vave 
Resistance of the equivalent system (kips) 

R Simplified notation of R - , the resistance function developed 

o A&--L 



between the g and (g-1) floors of a mult i -story building 

R The average value of R in the time interval from to t 1 
S g G. 

R . Resistance function developed between the g and the i floors 
gl of a mult i -story building (kips) 

R Maximum resistance of the g story of a multi-story building 
^ (kips) 

R Resistance developed in the spring which connects m to the ground 

S ^ & 

R , Total resistance acting on the g floor mass (kips) 

g"L 

R Maximum resistance developed by a structural system (kips) 

R . Plastic resistance of the "A" portion of two-way slabs where the 
edge is fixed (kips) 
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R 



'me 
*mf 
'n 
*u 

V 

*1 
3 

*lmA 
*lmB 



Plastic resistance of the "B !l portion of two-way slabs where the 
edge is fixed (kips) 

Maximum resistance of the equivalent system (kips) 

Fictitious maximum resistance (kips) 

Resistance of a structural element or system at time, t (kips) 

Horizontal ultimate load resistance of a shear wall (kips) 

Yield resistance of the structure (kips) 

Resistance in the elastic range (kips) 

Maximum resistance in the elastic range (kips) 

Maximum total resistance in the elastic range - "A" portion of 
two-way slabs (kips) 

MsLximum total resistance in the elastic range - "B" portion of 
two-way slabs (kips) 



R(t) Resistance as a function of time 
t 

g (t) ] 



H g (t) 



Resistance of the g story columns as a function of time 



R 



g 

s 

SE 

(SE) 

(SE) 



Time variation of the resistance of the equivalent single- 
5( ^ degree -of -freedom dynamic system for the g^ story 

Dosage of nuclear radiation without shielding (roentgens) 
Resistance as a function of displacement 

Radius of gyration of section (in.) 

Resistance per unit length of a beam or per unit area of slab 

Ratio of web reinforcement * A /bs 

Distance from point of explosion (ft) 

Roentgens 

Ratio of steel reinforcing placed perpendicular to the st.eel beam 

(in composite construction) in excess of that required to carry 

the slab bending stresses 

Average value of the resistance of the g story when the relative 
displacement between the g** 1 and the g-1** 1 story is negative 

Section modulus 

Strain energy absorption (kip -ft) 

Strain energy absorption of actual structure (kip -ft) 

Strain energy absorption of equivalent structure (kip-ft) 

Section modulus about the weak axis 

Coordinate apcis at k5 angle to x axis in x-y plane 

Spacing of stirrups and spacing of ties in reinforced concrete 

members (in*) 
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s Distance along s axis to centroid of loading (ft) 

s Distance along s axis to centroid of inertia force (ft) 

T Duration of the external load (sec) 
Resultant tension force (kips) 

T Fundamental period of vibration for complete spherical shell of 
c centroid thickness (sec) 

T Natural period of oscillation or fundamental period (sec) 

T , Lowest natural period of circular arch with pinned or fixed ends 

th 
T . Natural period of the i mode of oscillation (sec) 

T Time or rise of a step -load (sec) 



T-, T 9 -.., T Time duration of the external loads on the first, second, 
S +-u 

and g floors of a multi-story "building (sec) 

T 1X> T 2X*** T a:X Time d - uration when the effective loads on the first, 
^ second, and g* 11 stories are negative (sec) 

T 1 Period of oscillation of a fictitious system representing the 
811 g th story (sec) 

t Thickness of concrete slabs (in.) 
Thickness of flange of beams (in.) 

Time measured after the arrival of the incident shock wave (sec) 
Time variable 
Thickness of deep beam web 

t Time of arrival, or time required for the shock wave to travel 



a 



o^ir 
av 



from the point of explosion to the chosen location (sec) 
Average flange thickness of standard steel beam (in.) 



t. Time required for overpressure on the rear face of a closed 

rectangular structure to rise from zero to its maximum value (sec) 

t Time required to clear the front face of a structure from the 
reflection effects (sec) 

t. Time displacement factor; the time required for the shock front 
to travel from the frontmost element of a structure to the point 
or surface under consideration (sec) 

t e Time at which the limiting elastic deflection is reached 
t f Flange thickness of a WF steel beam (in.) 

t Time at which the limiting elastic deflection of the g story 
is reached 

t Time required for the vortex generated at the front face of a 
structure to travel a distance L f across the structure (sec) 
Time required for maximum displacement of element or structure 
to occur (sec) 

t Q Duration of the positive phase of the incident shock wave (sec) 
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t Time of rise; time required for the overpressure in the incident 
r shock wave to rise from zero to its maximum value (sec) 

t Thickness of stiff ener (in.) 
s 

t. Time required for the shock to pass over the length of the 
structure (sec) 

t Time of rise of an underground overpressure pulse (sec) 
t Web thickness of a steel beaon or channel (in.) 

t Time required to reach yield point of material (sec) 
yj? 

t... t ...t . t - Time sequence 
I 7 2 n rn-1 

t, , tp...t Time at which the maximum absolute displacements of the 
^ first, second and g^* 1 floor masses are reached 

tu* Time required for the overpressure on the exterior rear face of 
a rectangular structure with openings to rise from zero to its 
maximum value (sec) 

t 1 Time required to clear the front face of a structure with openings 
C from reflection effects (sec) 

t', t 1 Time at which the relative displacements X^ and X , respectively, 
are equal to zero 

t- Time immediately before t 
n n 

t+ Time immediately after t 
n n 

t. Time lag 
lag e 

At Time interval used in numerical analysis 

At * t ., - t A time interval 
n iw-1 n 

U Velocity of the incident shock front (fps) 

U. Velocity of the shock front of the shock wave formed in the 
interior of a rectangular structure with openings (fps) 

u Particle velocity in the incident shock wave (fps) 

Bond stress per unit of surface area of bar in reinforced con- 
crete design (psi) 

V Dynamic reaction (kips) 
Total shear (kips) 
Total vertical force on piles 

V. Total dynamic reaction along one edge "a," two-way slabs (kips) 

V Average vertical shear in length, L 1 (kips) 

V^ Total dynamic reaction along one edge u b/ ! two-way slabs (kips) 

V Total column load in flat slab design (kips) 

V Maximum shear capacity of web of deep beam section 

V Maximum vertical shear force (kips) 
max 
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v Vortex velocity (fps) 
Shear stress (psi) 

v Dynamic shear yield strength 

v = v (t \ Velocity of the g t floor at. time, t 
g,n g\ nj 

v Maximum shearing stress (psi) 
m* 

v = vft \ Velocity at time, t 
n ^ ny n 

v+ Velocity at t+ 
n n 

v- Velocity at t- 
n n 

v Horizontal velocity of axis of rotation "0" 
o 

vft ^ Velocity of mass at time t g 

v Static shear yield strength of steel (psi) 
^ Average unit shear stress in concrete (psi) 

v(t) Velocity as a function of time 

W Total energy yield of an atomic bomb expressed KT of TNT 
required for an equivalent total energy yield 
Weight (ibs) 

Total load on element of structure (kips) 
Work done (ft-lbs) 
Dynamic 1 , peak blast load 

W Work done on the actual system (ft-kips) 
a 

W Work done on the equivalent system (ft-kips) 

W Maximum work done on the g floor mass (ft-kips) 

W Absolute maximum work done on the g floor mass (ft-lbs) 
gm 

W m Maximum work done on the equivalent system by the equivalent load 

W Fictitious maximum work done on the equivalent system 

P 
W(t) Work done as a function of time 

w Uniformly distributed load (kips/ft) 

Length of channel shear connector (in.) 

Width of front face of rectangular building (ft) 

w n Weighting factor for vertical and batter piles 

X Relative displacement in a story of a multi-story building (ft) 

X Maximum relative displacement in the g story (ft) 

X Relative displacement of the mass g when the average external 
gs loads are applied statically 

X-,, Xp, X Relative displacement in the first, second, and g stories 
s of a multi-story building (ft), X x - x 

o o o"* 
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X Limiting elastic relative displacement for the g story column 
ge (ft) 

x Distance the incident shock wave travels after impinging on the 
frontmost element of a cylindrical segment or spherical dome to 
any point being considered 

Deflection of a structural element or system (ft) 
Distance of any pile from the centroidal axis 

x Collapse deflection 

x Limiting elastic deflection (ft) 

x Forced solution of a dynamic system for a given external load 

f'Vi 

x Absolute displacement of the g floor at time t 

g .k g 

x Maximum absolute displacement of the g floor (ft) 

S 01 th 

x -, x , x .. Absolute displacements of the g floor at time, 

g,n-J. g,n g,nn t t and t 

n-1' ir m-1 
x Maximum displacement (ft) 

Displacement at time t . , t , and t . , respectively 

tl""JL n IxT" _L 

x Initial displacement (ft) 

Horizontal displacement of axis of rotation "0" 

x Displacement of an elastic system subjected to the peak load B 
s acting statically (ft) 



4.) 



Displacement of a mass at time t 

4"V 

x Absolute displacement of the first, second, and g floors 
g of a multi -story building (ft) 

x Distance to centroid parallel to x axis (ft) 
x Horizontal velocity of axis of rotation "0" 

x Horizontal acceleration of axis of rotation "0" 

o ( 2\ 

x Acceleration in the x direction ^ft/sec j 

y Displacement or deflection of a structural system (ft) 

y Deflection of actual element (ft) 
a 

y Midspan deflection of actual element (ft) 
ac 

y Midspan deflection (ft) 

y Limiting elastic deflection (ft) 

Deflection of equivalent system (ft) 

y Limiting deflection in the elasto-plastic range (ft) 

y Maximum displacement (ft) 

y Displacement at time, t (ft) 
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y r Vertical distance from axis of rotation "0" to mass centroid of 
rotating mass, m 9 of structure 

y^ Deflection at midspan of column strip - flat slab (ft) 

y 2 Deflection at midspan of middle strip - flat slab (ft) 

y lm Maximum elastic range deflection (ft) 

y 2m Maximum elasto-plastic range deflection (ft) 

y Distance to the centroid parallel to y axis (ft) 

Vertical distance from axis of rotation "0" to mass centroid of 
total moving mass, m 

y Acceleration in the y direction (ft/sec J 

y Velocity in the y direction (ft/sec) 

y Velocity of displacement of the center of the beam (ft/sec) 

y Increment of displacement (ft) 

Z Plastic modulus of the cross section \in. / 

Z 1 Plastic modulus of the cross section in the weak direction 

z Depth of the resultant compressive force in concrete tee beam 

a Angle of incidence between the normal to the surface and the 
direction of propagation of the blast wave (degrees) 
Central angle of the arch 
Deflection coefficient for two-way slabs 
Design load ductility reduction factor 

a Angular acceleration of structure about axis of rotation l! !l 
(radians/sec J 

P Dimensionless ratio = 0.5^/1^ 7 
Ductility ratio 

6 Maximum bending deflection due to the application of a unit load 
b (ft/kip) 

& Shear wall lateral deflection at first cracking 
c 

5 Sub-area clearing factor for rectangular structures with openings 
n 

6 Maximum shear deflection due to the application of a unit load 

S (ft/kip') 

6 Shear wall lateral deflection at ultimate resistance 
u 

& Deflection at theoretical yield (ft) 

y ^ 

Angle of inclination to the horizontal of a gabled roof ( degrees j 
Parameter used to define the location of a point on the surface 
of a cylindrical segment or a spherical dome (degrees) 
End rotation of structural element 
Angular displacement 
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0_ Joint rotation at "bottom of column 

JD 

Rotation of beam at midspan (radians) 

G Angular displacement of structure about axis of rotation "0" 

6 m Joint rotation at top of column 

/ o \ 

H Viscosity of air in the incident shock wave \ Ib-sec/f t j 

Coefficient of friction 



roisson s ra^io 

Mass density of air in the incident shock wave I Ib-sec /in. ) 

/ 2/4 2, 4\ 
Mass density of soil I Ib-sec /in. , or Ib-sec /ft j 



v Poisson f s ratio 
P 

p s 

o Critical buckling stress (psi) 

$ Parameter used to define the location of a point on the surface 
of a spherical dome 

Phase angle in the transient response (radians) 
Internal friction angle (degrees) 

0. Phase angle of the i mode (radians) 
i 

7 Unit weight of soil (ibs/cu. ft) 

w Angular velocity of structure about axis of rotation "0" 
(radians/sec ) 
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ENGINEERING AND DESIGN 

DESIGN OF STRUCTURES TO RESIST THE EFFECTS OF ATOMIC WEAPONS 

WEAPONS EFFECTS DATA 
IMERODUCTIOIT 

3-01 PURPOSE AND SCOPE, This manual is one in a series issued for the 
guidance of engineers engaged in the design of permanent type military 
structures required to resist the effects of atomic weapons. It is appli- 
cable to all Corps of Engineers activities and installations responsible 
for the design of military construction* 

The material is based on the results of full scale atomic tests and 
analytical studies. The problem of designing structures to resist the ef- 
fects of atomic weapons is new and the methods of solution are still in the 
development stage. Continuing studies are in progress and supplemental 
material will be published as it is developed. 

The methods and procedures were developed through the collaboration 
of many consultants and specialists. Much of the basic analytical work was 
done by the engineering f irsn of Aramann and Whitney^ Hew York City, under 
contract with the Chief of Engineers. The Massachusetts Institute of Tech- 
nology was responsibley under another contract with the Chief of Engineers, 
for the compilation of material and for the further study and development 
of design methods and procedures. 

It is requested that any errors and deficiencies noted and any sug- 
gestions for improvement be transmitted to the Office of the Chief of Engi- 
neers, Department of the Army, Attention: ENGEB. 

3-02 REFERENCES. Manuals - Corps of Engineers - Engineering and Design, 
containing interrelated subject matter are listed as follows- 

DESIGN OF STRUCTURES TO RESIST THE EFFECTS 
OF ATOMIC WEAPONS 

EM 1110-3^5-^13 Weapons Effects Data 

EM 1110-3^5-^1^ Strength of Materials and Structural Elements 

EM U10-345-415 Principles of Dynamic Analysis and Design 



EM 1110-3^5-^13 3 ~ 02a 

1 July 59 

EM 1110-345-43*6 Structural Elements Subjected to Dynamic Loads 

EM 1110-345-417 SiJ3gle-Story Frame Buildings 

EM 1110-345-418 Multi-Story Frame Buildings 

EM 1HO-345-419 Shear Wall Structures 

EM 1110-345-420 Arches and Domes 

EM 1110-345-421 Buried and Semiburied Structures 

a. References to Material in Other Manuals of This Series, In the 
text of this manual references are made to paragraphs, figures, equations, 
and tables in the other manuals of this series in accordance with the 
number designations as they appear in these manuals. The first part of the 
designation which precedes either a dash, pr a decimal point, identifies a 
particular manual in the series as shown in the table following. 
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1110-345-^13 



1110-345-417 



paragraph 
3- 

5- 
6- 

7- 
8- 

9- 

10- 

11- 



flgure 

3- 
k. 

5- 

6. 

7- 

8. 

9. 

10. 

11. 



equation 


table 


it j 


3. 


(5. ) 


5- 


6. 




6. 


7- 




7- 


8. 




8. 


9. 




9- 


(10. ) 
11. ) 


10. 

11. 



ino-s^-^ 
1110-345-420 

1110-3^5-^21 

b Bibliography. A bibliography is given at the end of the text. 
Items in the bibliography are referenced in the text by numbers inclosed in 
brackets. 

c. List of Symbols. Definitions of the symbols used throughout this 
manual series are given in a list following the table of contents. 
3~3 RESCISSION* (Draft) EM 1110-3^5-^13 (Part XXIII - The Design of 
Structures to Resist the Effects of Atomic Weapons, Chapter 3 - Weapons 
Effects Data.) 

3-04 GENERAL. Atomic weapons may be exploded in the air, at ground level, 
or at various depths below the ground surface. The types of protected 
structures whose design are covered in this series of manuals are surface, 
shallow covered, and shallow buried structures. The major effects con- 
sidered in protection against atomic weapons are: 

1) Air blast 

2) Nuclear radiation 

3) Thermal radiation 

4) Crater ing 
(5) Ground shock. 

2 
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Underground burst and cratering cause insignificant effects on a 
structure unless the structure is in close proximity to the "burst. These 
effects have been given very limited coverage in this manual. 

The air burst is the major factor in the determination of the re- 
quired strength of a blast-resistant structure. To evaluate the loads on 
a structure it is necessary to develop the critical blast loading pattern, 
relative to each exposure of the surface of the structure. Pressure in- 
tensities will vary with respect to time, reckoned from the instant of 
arrival of the shock front. 

Nuclear radiation is an important consideration in the design of 
protective structures. In some designs radiation shielding may govern the 
required thickness of concrete and earth cover whereas the blast effects 
will be the determining factor relative to the designs for structural 
strength and stability. Consideration should be given to both direct and 
residual radiation. Incendiary properties and heat insulation from thermal 
radiation should be given attention in selecting suitable materials for 
construction of protective structures. 

3-05 BLAST WAVE FBENOMEMA. a. Infinite Atmosphere. To simplify the 
presentation of air blast phenomena, the explosion in an atmosphere of in- 
finite extent is first considered. Almost immediately after the detonation 
occurs, the expansion of the hot gases initiates a pressure wave in the 
surrounding air as represented roughly by the curve in figure 3-l(a). This 
shows the general nature of the variation of the overpressure (pressure 
above atmospheric) with distance from the explosion at a given instant 
during the early stages of shock formation. As the pressure wave moves 
outward from the center of the explosion, the following, or inner, part 
moves through a region which has been previously compressed and heated by 
the leading, or outer, parts of the wave. The pressure wave moves with the 
velocity of sound, and since this velocity increases with the temperature 
and pressure of the air through which the wave is moving, the inner part of 
the wave moves more rapidly and gradually overtakes the outer part, as 
shown in figure 3l(t>). The pressure wave front thus gets steeper and 
steeper, and within a very short period the pressure change becomes abrupt 
forming a shock front, as indicated in figure 3*l(c). At the shock front 
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Distance from Center of Explosion 
(o) 



Dittonct from Center of Explosion 
(b) 



..Shock 
Front 
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there Is a sudden increase of pressure from 
normal atmospheric to the peak shock pressure 
and the advancing shock front "becomes a moving 
wall of highly compressed air. 

Initially, in the hot central region of 
the explosion the pressure exceeds atmospheric 
"by a factor of many hundred thousand. The 
pressure distribution "behind the shock front 
is somewhat as illustrated in figure 3.l(c). 
It shows the peak overpressure at the shock 
front, dropping rapidly in a relatively small 
distance to a value about one -half the shock 
front overpressure. The magnitude of the 
overpressure is essentially uniform in the 
central region of the explosion. 

As the expansion proceeds, the pressure 
distribution In the region behind the shock 
front gradually changes. The overpressure is no longer constant but drops 
off continuously nearer the center. At later times when the shock front 
has progressed some distance from the center, a drop in pressure takes 
place below the initial atmospheric value at the center and a suction phase 
develops. The front of the blast wave weakens as it progresses outward 
and its velocity drops toward the velocity of sound in the undisturbed 
atmosphere. The sequence of events Just described is depicted in figure 
3-2. This shows the overpressure dis- 
tribution in the blast wave as a func- 
tion of the distance from the explosion 
at different stages in the expansion. 
When the suction phase is well de- 
veloped, the overpressure in the blast 
wave resembles the heavily drawn curve 
in figure 3.2. The pressure variations 
in the blast wave from this time on can 
be considered in two different ways as 



I 



Oistonct from Center of Explosion 
(c) 

Figure 3. L Overpressure during 

early stages of shock 

front formation 




Distance from Center of Explosion 



Figure 3.2. Variation of overpressure with, 
distance at various times 
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shown in figures 3-3 and 3Aa. 

First, the heavy curve of figure 
3.2 has been redrawn in figure 3-3 to 
show the variation of overpressure 
with distance at a given time. The 

symbol P represents the peak 
so 

overpressure, or shock intensity, in 
pounds per square inch and U repre- 
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Dittonce from Center of Explosion 

Figure 3.3. Variation of overpressure with 
distance at a given time 



sents the velocity of the shock front 
in feet per second. The arrows under the curve show the direction of move- 
ment of the air mass, or blast wind, in the positive and negative phases. 

The peak overpressure in the positive 
phase is higher than the maximum 
overpressure in the negative phase. 
Consequently, the blast wind is of 
higher velocity in the positive phase 
than in the negative, or suction, 
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Time after Detonation 



Figure 3.4a. Variation of overpressure with 

time at a given location phase . 

Second, the pressure variations in the same wave may be considered 
at a point located a given distance from the explosion as indicated in 
figure 3*^-a, which shows the variation of overpressure with time at a fixed 

location. The symbol t is the time of arrival, or the time in seconds 

a 

for the shock front to travel from the center of the explosion to the given 
location, t is the duration in seconds of the positive phase, and p 
is as previously defined. 

The rate of decay of the positive pliase varies with respect to the 
magnitude of the peak overpressure* In figure 3-^b the rate of decays 
for various magnitudes of peak overpressures is expressed in terms of 

the ratio P/P to the ratio t/t , where P is the overpressure at 

7 so ' o* s * 

time t after the arrival of the shock front. The peak negative over- 
pressure is approximately one-eighth of the peak positive overpressure 
and the duration of the negative phase is approximately four times the 
duration of the positive phase of % the blast. 

The peak overpressure P in free air before the shock front 

so 
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Figure 3.4b. Overpressure decay 
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Figure 3.4d. Dynamic pressure decay 
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Slant Distance from Burst (ft) 

Figure 3.5. Peak overpressure for 1-KT free air burst by 
2W assumption (taken from reference[36]) 
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Center of 
Explosion 



Radius of 
Shock Front 



strikes the ground is a function of the distance from the point of "burst 
and the yield of the weapon. Values for 1-KT bomb given in figure 3-5 are 
from curves given in reference [36], which are based on the supposition that 
a contact surface burst of W kilotons energy is equivalent in blast 
characteristics to an explosion of 2W kilotons high in the air and prior 
to reflection. 

b. Finite Height of Burst* The discussion in the preceding para- 
graph deals with the air blast from an atomic bomb exploded in an atmos- 
phere of infinite extent. A discussion follows relative to the influence 
of the surface of the earth on propagation and attenuation of the air blast 

If the bomb is detonated at a height h above the surface of the 
earth, the shock front will have the general configuration indicated in 

figure 3.6 during the brief time inter- 
val before it impinges upon the surface., 
The shock front is the same as if propa- 
gated in an atmosphere of infinite ex- 
tent. 

A short time later the radius of 
the shock front becomes greater than 

r . - Ol . f , . ... h, and a portion of the shock front 
figure 3.6. shock front before sinking 

ground surface impinges upon the earth ! s surface and 

is reflected back, forming the reflected shock front illustrated in 
figure 3-7* The arrows at the incident and reflected shock fronts indicate 
the direction in which the blast waves 
are traveling. The location of the re- 
flected shock front is roughly deter- 
mined by drawing an arc with center lo- 
cated a distance h below the earth's 
surface and directly under the point of 
detonation, joining the points of inter- 
section of the incident shock front with 
the earth's surface. The symbol a 




/ I \ Reflected 
a . a \ Shock 
Front 

Earths 
^Surface 




Figure 3. 7. Shock front reflection at 
ground surface, a less than 45 



represents the angle of incidence of the shock front with the earth 's sur- 



face. The reflected shock front overpressure P 



j- \ 



is a function of the 
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incident shock overpressure P and the angle of incidence a [20]. 

so 

The reflected shock front in figure 3-7 travels through the atmos- 
phere at a higher velocity than the incident shock "because it is traveling 
through a region of greater than atmospheric pressure and gradually over- 
takes and merges with it to form a single shock front called the Mach Stem 
as shown in figure 3*8. The fused shock front thus formed is normal to and 
travels parallel to the earth's 
surface. The Mach Stem forma- 
tion is initiated when the 
angle of incidence a of 
the shock front "becomes greater 
than approximately 45 degrees. 
Once formed, the height of the 
Mach Stem gradually increases 
as the radius of the shock 
front "becomes greater. 

The region on the earth's 
surface within which a is 
less than 45 degrees, and no 
Mach Stem is present, is called 



Region of 
Regular Reflection 



Region of 
Mach Reflection 




Figure 3.8. Shock front reflection when & is 
greater than 45 



the region of regular reflection, while the region for which a is greater 
than 45 degrees and a Mach Stem is present is called the region of Mach re- 
flection. The importance of the Mach Stem phenomenon is that it results 
from the fusing of two shock fronts to form a single shock front of higher 
overpressure and of greater destructive potential to structures located in 
its path. 

The shock front velocity U is a function of the peak overpressure 

P ; its value for standard atmospheric conditions is plotted in figure 3-9- 
so 

The time of arrival t of the shock front at a given location is a func- 

a 

tion of the distance from the center of the explosion and the total energy 
yield of the weapon. The duration on the ground surface of the positive 

phase t of the blast wave is a function of the peak overpressure P 

o so 

and the total energy yield of the weapon. Values of the duration of the 
positive phase are plotted in figure 3.10a and in figure 3.10"b. 
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Figure 3.10a. Duration of positive phase typical air burst 
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Figure 3.10b. Duration of positive phase surface burst 
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Reflected overpressure ratio P /P is plotted in figure 3.11 as 
a function of the angle of incidence a of the shock front. This figure 
applies to an inclined shock front striking the earth's surface or to a 
vertical shock front striking a reflecting surface such as a wall of 
structure . 

The peak overpressure P in the "blast wave (shock front overpres- 

so 

sure) at the ground surface plotted against distance from ground zero for 
various bomb sizes is given in figure 312a for a typical air burst and in 
figure 312b for a surface burst* The typical air burst is at a height 
above the ground which will maximize the target area covered by the 12-psi 
blast pressure level* This height probably cannot be attained by manned 
aircraft for the larger size bombs * 

For other heights of burst P may be obtained from charts in 

s o 



reference ^9 ? where the peak overpressure is plotted as a function of the 
height of burst and the horizontal distance from ground zero* 

The variation of the overpressure in the blast wave on the ground 
with time can be determined in the same manner as for the bomb burst in an 
infinite homogeneous atmosphere* Equation 3-1 plotted as figure 3^d is 
valid for ideal wave shapes which occur at approximately 10 psi or lower* 

-t/t 

P = P (1 - t/t )e (3*1) 

s so ' o 

For higher overpressure levels the variation of overpressure with time at 
the ground level is usually irregular^ due to thermal effects and to other 
modifying influences of the ground surf ace . Reference [^9] contains data 
on these irregular wave forms. However, due to lack of data and the com- 
plexity of calculations involving irregular wave forms^ the ideal wave forms 
that are given in figure 34c should be used for practical design purposes. 

The dynamic pressure q due to the motion of the air particles in 
the blast wave is given by the relation q = pu^/2 where p is the mass 
per unit volume of the air and u is the velocity of the air particles. 
The peak values of q at the ground surface plotted against distance from 
ground zero for various bomb sizes are given in figure 323b The peak 

values of q plotted against the peak overpressure P are given' in 

so 

figure 3- 23 a 
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The rate of decay of the dynamic pressure is a function of the peak 
overpressure* Figure 34d shows variation in dynamic pressures for ranges 

of overpressures as a function of q/q, at a fraction of time t/t ; 

so o 

where P is the over-pressure and CL is the peak dynamic pressure and 
so o 

q_ is the dynamic pressure at time t , after the arrival of the shock 
front . 

3-06 SCALING BLAST PHENOMENA. It has been found that air blast phenomena 
such as the pressure and duration at different distances are related, for 
different strength bombs , according to the ratio of the cube root of the 
yield expressed in TNT. These relations are referred to as scaling laws. 

These scaling laws state that if a given peak overpressure is experi- 
enced at distance r, from an explosion of a bomb of total energy yield 
W- the same peak overpressure will be experienced at distance r from 

lelc 
1/3 

(3-?) 



the explosion of a bomb of total energy yield W p where: 




The same scaling laws also state that while the peak pressures from 
the two bombs are equal to the two radii r.. and r p , the durations of 
the blast pressure waves at the two points are different. If the duration 
of the positive phase of the pressure wave from the first bomb is t . at 
radius r.. the duration of the positive phase t of the pressure wave 
from the second bomb at distance r p will be: 




LOADING ON EECTANGUIAR STRUCTURES WITHOUT OPENINGS 

3-0? LOADING ON STRUCTURES. The manner in which the blast wave loads a 
structure varies with the ratio of the distance of the structure from 
ground zero to the height of the burst. The loading on a structure located 
within the region of regular reflection (figure 3.13) is greatly different 
in character from that for a structure located beyond this region. Since 
little or no data are available within the region of regular reflection and 
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Reflected 
Shock Front 



Incident 
Shock Front 
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Shock Front 




Structure 
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Figure 3. 13. Structure located in region 
of regular reflection 



Figure 3. 14. Structure located in region 
of Mack reflection 

since a surface structure located in this region would be subject to such e 

high overpressure as to preclude a practical design, the procedures for th* 
determination of loads are restricted to surface structures located outside 
the region of regular reflection (figure 3.14) where the Mach Stem is high. 
enough to cover the structure 

The problem is that of computing the loading on a structure due to 
the impingement upon it of a blast wave traveling parallel to the surface 
of the earth* The structure is considered as being oriented with one face 
normal to the direction of propagation of the shock wave^ since such an 
orientation produces the most severe loading on the structural elements* 
For the design of a structure so located^ it is necessary to predict the 
overpressures which would exist on various portions of the structure as a 
function of time t measured after the shock front strikes the front wall 

The procedures developed in the following paragraphs for the deter- 
mination of the loads on different types of structures are based on exper 
mental results obtained from TM 1 and atomic explosions and from shock tube 
data. The computational methods presented are based on an interpretation 
of these experimental results and the equations given in connection with 
the loadings on a structure result from curve- fitting of measured data and. 
should not be regarded as representing the development of any new theory , 
Nevertheless^ every attempt has been made to make the empirical relation- 
ships consistent with such theory as now exists. 

The following paragraphs are devoted to a general presentation of tli 
procedures for the determination of the loads on structures of various 
types , Pertinent quantities are introduced and explained^ and justificati 
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for the methods employed is given as well as the bases for their concep- 
tion. The general presentation for each type of structure is followed "by a 
paragraph giving a step-by-step procedure by which the loads can be com- 
puted. Illustrative numerical examples are presented for the closed 
rectangular structures* 

In using procedures presented herein^ inconsistencies in loads com- 
puted may be encountered^ particularly in the case of buried and semiburied 
structures and with peak overpressures greater than 25 psi* The lack of 
experimental data and an adequate theory in these regions requires use of 
engineering judgment and some degree of conservatism* Extrapolation and 
use of load data for peak overpressures greater than, 25 psi are strictly 
an expediency and must be recognized as such by the engineer* 

a. Diffraction and Drag Loading^ The loading on aboveground struc- 
tures resulting from the air blast produced by an air or surface burst may 
be considered to consist of a diffraction phase and a drag phase . 

The diffraction phase of the loading is the term given to the initial 
phase of the blast loading on a structure when the reflected pressures as- 
sociated with the air blast are acting on the structure. The time required 
for the blast wave to surround the structure completely and the presence of 
large reflected pressures on the front wall cause net lateral loads to be 
exerted on the structure as a whole in the direction of travel of the blast 
wave. The local and differential forces which act on the structure during 
the initial stages are defined as the diffraction phase loading. 

The drag phase of the loading is the term given to the second phase 
of the loading on a structure due to the mass and velocity of the air 
particles in the blast wave after the envelopment of the structure by the 
wave front and the reflection effects have decayed. This phase of the 
loading is most important on open structural frameworks and on structures 
having small dimensions such as stacks,? poles,, chimneys^ etc* 

In general, the diffraction phase of the loading can be neglected and 
only the drag phase considered if the minimum dimension of the structure 
perpendicular to the direction of travel of the shock wave and the dimen- 
sion of the structure parallel to the direction of travel of the blast 
wave are less than approximately 5 ft. Because of such short dimensions^ 
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the duration of the diffraction phase is very short* The larger the loaded 
area the greater is the time required for the transient reflection effects 
to decay and reach the relatively steady condition of the drag phase * The 
longer the duration of the reflection effects the more important is the 
diffraction phase of the loading . In all cases except when a precursor 
exists,, the pressures on any portion of the structure existing during the 
diffraction phase are greater than during the drag phase * 

Large quantities of heated dust may be raised by the thermal radia- 
tion before the arrival of the blast wave* The air near the ground is als- 
heated by the thermal radiation and this may cause a precursor wave to fon 
ahead of the main shock* Test results indicate that when precursors form, 
there is usually a dust-laden atmosphere present and the dynamic pressures 
are much larger than when there is no precursor. These high dynamic pres- 
sures are usually observed in the high overpressure regions. 
3-08 LOADING ON CLOSED RECTANGULAR STRUCTURES. The type of structure for 
which the loading predictions of this section are applicable is illustrate 
in figure 3*15* The behavior of the blast wave upon striking a closed rec 
tangular structure is depicted in figure 3l6(a), (b)^ (c), and (d). This 
figure illustrates the position of the shock front and the behavior of the 
reflected and diffracted waves at successive times during the passage of 
the blast wave over the center, portion of the structure of figure 3*15- & 
the shock front strikes the wall of the building (figure 3*l6(a)) a re- 
flected blast wave is formed and 
the overpressure on this wall is 



-*- Incident 
Shock Front 



raised to a value in excess of the 
peak overpressure in the incident 
blast wave* This increased 
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Figure 3. 15. Closed rectangular 
structure 



Figure 3* 16. Behavior of blast wave along cen- 
ter portion of closed rectangular structure 
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overpressure is called the reflected overpressure and is a function of the 
peak overpressure in the incident blast wave and the angle of incidence of 
the shock front with the front wall which is zero degree in this case. At 
the instant the reflected shock front is formed, the lower overpressure 
existing in the incident "blast wave adjacent to the top edge of the front 
wall initiates a rarefaction wave (figure 3l6("b)), or a wave of lower 
ove rp res sure than that which exists in the reflected "blast wave* This 
rarefaction wave travels with the speed of sound in the reflected "blast 
wave toward the "bottom of the front wall. Within a short time, called the 
clearing time, the rarefaction wave enfeebles the reflected blast wave to 
such an extent that its effect is no longer felt and reduces the overpres- 
sure existing on the front wall to a value that is in equilibrium with the 
high-velocity air stream associated with the incident blast wave. The 
overpressure on the front wall when equilibrium with the high- velocity air 
stream is reached is equal to the stagnation overpressure at the base of 
the front wall and an overpressure somewhat less than that in the blast 
wave at the top edge of the front wall. The stagnation overpressure is de- 
fined as that overpressure existing in a region in which the moving air has 
been brought completely to rest causing the pressure intensity to be in- 
creased by the loss in momentum of the air particles. 

At some time after the shock front strikes the front wall of the 
structure, equal to the length of the structure divided by the shock front 
velocity, the shock front reaches the back edge of the structure and starts 
spilling down toward the bottom of the back wall (figure 3*l6(c)). The 
back wall begins to experience increased pressures as soon as the shock 
front has passed beyond it. The effect is first observed at the top por- 4 
tions of the back wall and proceeds toward the bottom. A vortex, which is 
a region of air spinning about an axis at a high speed with low overpres- 
sures existing at its center because of the venturi effect, is created on 
the back wall and grows in size, traveling toward the base from the top 
edge and also moving away from the wall (figure 3l6(b) and figure 317)* 
The maximum back wall overpressure develops slowly as a result of the vor- 
tex phenomena and the time required for the back wall to be enveloped by 
the blast wave. 
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Figure 3.17. Photograph of a blast wave passing over a square block. 
Light and dark bands represent contours of constant pressure 

As the shock front passes beyond the front vail (figure 3.l6(b)) the 
overpressure exerted on the roof of the structure is initially raised to a, 
value nearly equal to the overpressure existing in- the incident "blast wav . 
However, the air flow caused by the pressure difference between the re- 
flected overpressure on the front wall and the blast wave overpressure on 
the roof causes the formation of a vortex along the top edge of the front 
wall. The vortex travels with a gradually decreasing intensity along the 
roof of the structure (figure 3*l6(c)) at a slower rate than the shock 
front velocity. It causes a decay of the overpressures built up by the in- 
cident blast wave. After the passage of this vortex, the higher overpres- 
sures in the blast wave again become dominant and cause a second build-up 
of overpressures along the roof. 

If a horizontal section through the structure is considered as shown 
in figure 3l8> it will be found that the blast wave produces effects simi- 
lar to those produced on the roof, front, and back walls as indicated in 
figure 3*l6- The major effects which should be noted are: 
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(1) The formation of rarefac- 
tion waves at each end of the front 
wall of the structure which travel 
toward the center to relieve the 
higher reflected pressures on the 
front wall; 

(2) The diffracted wave 
fronts at the rear edges of the side 
walls which travel toward the center 
of the back wall; 

(3) The vortices which form on 
the side walls near the front edges 
and travel slowly toward the "back; 

(4) The vortices which form at 
the ends of the back wall and travel 
slowly toward its center* 

In determining the loads on a 
structure, it is convenient to use 
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Figure 3.18. Behavior of blast wave along horizon-" 

,, . , , , , . , , n T _ , iaZ section through closed rectangular structure 
the instant at which the shock front 6 6 

impinges on the front face as a reference for time (t =0). In adopting 
this convention, it is necessary to introduce a time -displacement factor 
t , which is the time required for the shock front to travel from the front 
face of the structure to the surface or point under consideration. Figure 
3.19 illustrates this convention. In this figure, t - t , is the time 
after the shock front has passed a given location on the structure and de- 



Time Scale 




Figure 3. 19. Time-displacement 
factor convention 



termiries the overpressure P in the 

s 

incident blast wave at that location. 
For example,, the time-displacement fac- 
tor for a given location on the roof a 
distance L T from the front face of 

the structure is t, = L T /U 

a ' o 

The total loading on any face of 
a cubical structure is equal to the 
algebraic sum of the respective 
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overpressures or reflected overpressures combined with the dynamic pres- 
sures. The later pressures are due to the wind effects of suspended parti- 
cles and are referred to as If drag M loading which may result in either posi- 
tive or negative corrections depending on the orientation of the surface 
with respect to the direction of travel of the shock front. 

a. Average Front Wall Overpressure P .At the moment the inci- 

i ron ~c 

dent shock front strikes the front wall, the overpressure on the front wall 
is immediately raised from zero to the reflected overpressure P -, Re- 
flected overpressure is plotted in figure 320 for zero angle of incidence 
as a function of the peak overpressure of the incident shock front. The 

value of P -, from figure 3.20 is the same as P for a = using 

rex j- j.*"*uc 

figure 3.11. The incident shock front continues its motion over the top oJ 
the structure, while the reflected shock front moves away from the front ol 
the building in the opposite direction. Initially, the air pressure be- 
tween the reflected shock front and the front wall is higher than the pres- 
sure "behind the incident shock front. This causes air to move around to 
the sides and over the top of the structure into the lower pressure zone "b 
hind the incident shock. The rarefaction waves thus formed move from the 
end and top edges toward the center of the front face with the speed of 
sound for the pressure existing in this region of reflected overpressure [l 
For the range of shock strengths used in the Princeton shock tube ex 
periments [^] which correspond to peak shock overpressures of 2 to 50 psi 
in a standard atmosphere, it has been found that the time required to clea 
the front face of reflection effects is determined by the dimensions of th 
front face and the peak overpressure of the incident shock wave. This 
clearing time t is given by the relation: 

where 

h r = clearing height, taken as the full height of the front face 
or half its -length whichever is the smaller 

refl = velocit y o:f sound in the reflected region, plotted as a func- 
tion of the peak overpressure of the incident shock wave in 
figure 3-21 

During the time required to clear the front wall of reflection ef- 
fects, the average overpressure on the front wall decreases from the 

26 



3-08a 



EM 1110-345-413 
1 July 59 



320 




20 40 60 80 100 

Peak Incident Overpressure, f| (psi) 

Figure 3.20. Reflected overpressure vs peak incident overpressure for normal reflection 
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reflected overpressure to a value given "by the following equation: 
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front 



= P + 0.85<i 
s 



(3.5) 



where 



P = overpressure in the incident "blast wave at the front vail at 
any time t - t , as given by the equation 



P /P 
7 



s so 



-, -(t - 1 )/t 

-<t-t a )/t ]. d 



(3-6) 



Equation (3-6) is plotted in figure 3.22, and tabulated in table 3.1, where 

the ratio P /P is given as a function of (t - t,)/t This equation is 
s so d o 

similar to equation (3-1) with (t - t,) substituted for t. 

In paragraph 3-05"b the dynamic pressures associated with high over- 
pressures were discussed. These pressures were evaluated in figure 323a. 
The rates of overpressure decay and dynamic pressure decay were set forth 
in figures 3.4b and 3.4c. 

Design requirements for aboveground frame structures are usually 
for peak overpressures of 10 psi or less* The use of equations (3*6), 
(3*7a), and 3.7c following and tables 3*1 and 32 will be applicable 
and convenient* 

The quantity q in equation (3-5) is the dynamic pressure due to the 
motion of air particles in the incident blast wave at any time (t - t,) and 
is given by q = pu^/2 where P is the mass per unit volume of the air and u 
is the velocity of the air particles. Values of dynamic pressures vs over- 
pressures for the theoretical (without precursor) and precursor condition 
are given in figure 3-23a. It will be noted that the values are about the 
same for each case up to about 10-psi overpressure. The dynamic pressure 
vs distance for miles for various weapon sizes is given in figure 323"b. 



q. = 



Peak dynamic pressure q 



+ i (P B /lA.7)_ 
for use when overpressures are less than 



(3-7a) 



10 psi or where there are no precursor effects considered is given by the 
equation 

'Z i " \ Jr / JL.T" I I 

(3.7^) 



1 (P an /lA.7) 



so' 
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Table 3. L Blast Wave Overpressure Ratio vs Time Ratio 



(t - t d )A 


P /P 
s' so 


(t - t d )/t o 


P /P 
s' so 


(t - t d )/t 


P /P 
s' so 


0.00 


1.000 


0.30 


0.519 


0.60 


0.220 


0.01 


0.980 


0.31 


0.506 


0.6l 


0.212 


0.02 


0.961 


0.32 
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0.62 


0.204 


0.03 


0.941 


0.33 


0.482 


0.63 


0.197 


0.04 


0.922 


0.34 


0.470 


0.64 
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0.05 


0.904 


0.35 


0.458 


0.65 


0.183 


0.06 


0.885 


0.36 


0.447 


0.66 


0.176 


0.07 


0.867 


0-37 


0.435 


0.67 
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0.08 


0.849 


0.38 


0.424 


0.68 
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0.09 


0.832 


0.39 
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0.69 
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0.10 
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0.40 


0.402 
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0..149 


0.11 
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0.4l 
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0.71 


0.143 


0.12 
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0.42 
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0.72 
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0.13 


0.764 


0.43 


0.371 


0.73 
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0.14 
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0.44 


0.361 


0.74 


0.124 


0.15 
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0.45 
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0,75 
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0.16 


0.716 


0.46 
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0.76 


0.112 


0.17 
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0.22 
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0.612 


0.53 


0.277 


0.83 


0.074 


0.24 


0.598 


0.54 


0.268 


0.84 


0.069 


0.25 
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0.55 
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0.064 


0.26 
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0.56 


0.251 
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0.27 


0.557 


0.57 
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0.87 
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0.28 


0.544 


0.58 
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0.88 


0.050 


0.29 


0.531 


0.59 


0.227 


0.89 
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Figure 3.23a. Peak dynamic pressure vs peak incident overpressure 
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To obtain g as a function of time for overpressures less than M 
10 psi it is convenient to use the ratio qJq. Q , plotted in figure 3- 2k and 
tabulated in ta"ble 3.2 as a function of (t - "t d )/t Q , from the following 



expression 



-3-5(t - t d )/t c 



(3-7c) 




0.2 0,4 0.6 0.8 

Time Ratio,(t- t d ) /t 



I.O 



Figure 3.24. Dynamic pressure ratio vs time ratio 
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<* - tcA 


*/<i 


(t - t d )t 


*/*o 


<* - V^o 


<AO 


0.00 


1.000 


0.30 


0.245 


0.60 


0.049 


0.01 


0.956 


0.31 


0.233 


0.61 


0.046 


0.02 


0.914 


0.32 


0.222 


0.62 


0.043 


0.03 


0.873 


0.33 


0.211 


0.63- 


0.041 


0.04 


0.835 


0.34 


0.201 


0.64 


0.038 


0.05 


0.798 


0.35 


0.191 


0.65 


0.036 


0.06 


0.762 


0.36 


0.182 


0.66 


0.034 


0.07 


0.728 


0.37 


0.173 


0.67 


0.032 


0.08 


0.695 


0.38 


0.164 


0.68 


0.030 
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0.664 


0.39 


. 156 
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0.028 
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0.70 


0.026 


o.n 
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0.4l 


O.l4l 
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0.024 
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0.42 
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0.72 
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0.13 


0.552 


0.43 


0.127 


0.73 


0.021 


0.14 


0.527 


0.44 


0.120 


0.74 


0.020 


0.15 


0.503 


0.45 


o.n4 


0.75 


0.018 


0.16 


0.480 


0.46 


0,108 


0.76 


0.017 


0.17 


0.458 


0.47 


0.102 


0.77 


0.016 


0.18 


0.437 


0.48 


0.097 


0.78 


0.014 


0.19 


0.417 


0.49 


0.092 


0.79 


0.013 


0.20 


0.397 


0.50 


0.087 


0.80 


0.012 


0.21 


0.379 


0.51 


0.082 


0.81 


0.011 


0.22 


0.361 


0.52 


0.078 


0.82 


0.010 


0.23 


0.344 


0.53 


0.074 


0.83 


0.009 


0.24 


0.328 


0.54 


0.070 


0.84 


0.008 


0.25 


0.313 


0.55 


0.066 


0.85 


0.008 


0.26 


0.298 


0.56 


0.062 


0.86 . 


0.007 


0.27 


0.284 


0.57 


0.059 


0.87 


0.006 


0.28 
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0.58 
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0.006 
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t - t. 



where 



For the front vail overpressures^ t, = 0, hence w v 
t measures the time after the instant at which the shock front impinges on 
the front wall. 

Using the above quantities^ the curve of average front wall overpres- 
sure P , versus time may be determined. Figure 3-25 shows a typical 
front 

front wall average overpressure curve. The curve a-b-c is defined by 

equation (3-5). Point d, the re- 
flected overpressure^ is connected by 
a straight line to point b, the aver- 



r rtfl 



Note; 

<s*P so +0.85q 

= 3h'/c r @fi 




3/4t c t c 3/2t c 



Time, t 

Figure 3.25. Average front wall overpressure 
vs time*" closed rectangular structure 



age overpressure at time t - t , * t . 
Since the resulting discontinuity at 
point b is not compatible with 
actual behavior of the loads on the 
front face^ these two curves are 
smoothed by fairing in curve e-f as 
shown. The curve of average over- 



pressure versus time on the front wall is then defined as curve o-d-e-f-c. 

The foregoing method for determining the curve of average overpres- 
sure versus time on the front wall is based on the assumption that the time 
required for the overpressure in the incident shock wave to rise from zero 
is negligible. 

If this is not the case^ the time required for the peak overpressure 
to rise from zero^ called the time 
of rise t , must be taken into 
account. This is accomplished by 
first computing the average front 
face overpressure for a zero time 
of rise as described above * Curve 
o-d-e-f-c of figure 3,26 illus,- 
trates the average front -wall 
overpressure computed for a zero 
time of rise. Locate point g 

with time coordinate t and over- 

r 




Tlme, t 



Figure 3.26. Graphical procedure for evaluating 

reduced average front wall overpressure due 

to time of rise of incident shock wave 



pressure coordinate P 



refl* 



Join points o and g vith straight line 
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intersecting d-e-f at h. The average front wall overpressure for time of 
rise t is given by curve o-h-e-f-c. 

The overpressure on the front vail of a structure is not uniformly 
distributed. The maximum value occurs at the mid-point of the base and the 
minimum value occurs along the edges* Those portions of the front face 
nearest to the edges are cleared of the reflection effects in a shorter 
time than the remainder of the front face and the overpressure existing at 
those points is lower, following the clearing stage. The net effect of 
this vertical and horizontal variation is negligible and of questionable 
value for design purposes . Hence, the front vail loading is assumed to be 
distributed uniformly over the front wall surface. For the same reasons, 
the rear vail loading given next is assumed to be uniformly distributed 
over the rear vail. 

b. Average Back Wall Overpressure P , . For the back wall, the 
time-displacement factor t is ' L/U , where L is the length of the 
building in the direction of propagation of the shock and U is the 
shock front velocity* 

When the shock front crosses the rear edge of the structure, the foot 
of the shock spills down the back vail. The overpressures on the back vail 
behind this diffracted vave are considerably less than those in the inci- 
dent blast vave due to the vortex vhich develops at the top and travels 
down the vail- A period of time longer than that required for the travel 
of this diffr&cted shock to the bottom of the vail must pass before the 
back vail average overpressure reaches its peak value The time required 
for this build-up to occur, t, , measured from the instant at vhich the 
Shockwave reaches the back wall is equal to ^h f / c 9 vhere h f is the 
clearing height of the back wall, taken equal to either the full height 
of the back vail or half the width of the building, whichever is the 
smaller, and c is the velocity of sound in undisturbed air, 1115 fps. 

The peak value of the average overpressure on the back wan after 
this build-up has been completed is : 



' P sb <V2) l + (1 - *- (3.8) 

where P , is the incident blast vave overpressure at the back wall at 

time t - t, t, j (P, . ) is the peak value of the average overpressure 
d D' v back'max * * 
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t. + t ; 3 

0. D 



0.5P 



SO 



on the back wall which occurs at time t 
e SB 2.7183 = base of natural logarithms. It is assumed that P gQ in the 
incident blast wave does not diminish in strength as the wave passes over 
the structure. 



Figure 3.27a illustrates the variation of the ratio P, 



bac 



k /P s with 

time. Per times in excess of t * t + t, the relation is given by the 
following equation: 



\)ack 


(P ) 

back max 


1 - 


(p ) 

v back'max 


ft - (t d + V 


? 3 ' 


P st + 


P sb 


L t o~ \ 



(3.9) 



where t is the duration of the positive phase. 



Figure 3.27b shows a typical back wall average overpressure curve, 




V L/U 

t b =4h/c 



/ Time, t 

For t d < t <t b -f t d , P |XJck Vories Lineorly 
from Zero at t=t d .to the Value Indicated 
ot tl b +t d . 

Figure 3.2? 'a. Average back wall 
overpressure ratio vs time 
closed rectangular structure 



Note: Fnet"Pfront-Pback 

Average Front Wall 
Overpressure, Pf ron t 



Average Back Wall 
Overpressure, P bach 



Net Horizontal 
Overpressure, P net 





Time, t 

Figure 3.27b. Average back wall overpressure 
vs timeclosed rectangular structure 

c. Average Net Horizontal Over- 
pressure P_^^_. Considering as positive 



all overpressures exerted on the struc- 
ture and directed toward the interior, 



net 



front 



back 



(3-10) 



where 



are 



Time, t 



Figure 3.28. Net horizontal overpressure 
timeclosed rectangular structure 



vs 



net' front' w ,~xv 
given in terms of time t. Figure 3-28 

illustrates graphically the variation 
of the net average horizontal over- 
pressure with time. 
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d. Local Roof Overpressure P . The procedure developed for pre- 
diction of overpressures on the roof of a structure is based primarily on 
curve-fitting of the Princeton and Michigan shock tube data [4 through 13]. 
Pressures calculated by this procedure have been found to agree satisfacto- 
rily with overpressure records obtained in the GREENHOUSE structure testing 
program and the Sandia HE tests [6j. 

During the passage of the blast wave across the structure, low pres- 
sure areas develop on the roof and side walls due to vortex formation as 
indicated in figure 3.16 and figure 3*18. Shock tube data indicate that 
this vortex detaches itself from the front edges and moves across the 
structure with gradually increasing 
speed. Vortices are formed all along the 
front edges of the roof and sides of the 
structure, that is along a-b-c-d of 



Area Affected 
by Presence 
of Vortices 




figure 329. At corners b and c 

where they are aligned at 90 to each 

other, the vortices tend to interfere 

with each other and to move away from 

the roof and wall surfaces. This results in a diminution of the effects 



Figure 3.29. Areas on roof and sides 
most affected by vortex action 



of the vortices at these edges which proceeds along the axes of the vor- 
tices from corners b and c as the vortices move toward the rear of the 
roof. Consequently the regions on the roof and side walls which are 
strongly affected by the vortices do not extend to all edges of the roof, 
but are triangular or trapezoidal in shape as indicated in figure 3-29. 




C.L. 



y^- Length, L 



These observed results indi- 
cate a variation of overpres stares on 
the roof in a direction parallel to 
the shock front in addition to a 
variation along the roof in the di- 
rection of propagation of the wave. 
Although this lateral variation is 
a smooth one, the roof and side 
walls have been divided into three 



Figure 3.30. Location of loading zones on 
roof and sides of structure 

zones as illustrated in figure 3-30 to facilitate computation of local 
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pressures . Although only the location of the zones on the roof is dis- 
cussed below ^ the zones on the side wall are similarly located as shown in 

figure 330 

Zone 1. This zone is a strip on the roof extending from the sides 

toward the center line a distance equal to one -quarter of the length L of 



the building. 



Blast Wove 



->i C.L. 



1.0 



o 
K 




Note: t d L'/U 



Time, t 



t d *t 



Figure 3.31. Local roof or side wall over*- 

pressure ratio vs time. Zone Inclosed 

rectangular structure 



for any point in Zone 1 is equal to 
the overpressure-time curve for the 
incident blast wave displaced in 



time by factor t , 



This 



d i_j i i^ o 
o 

is consistent with observed vortex ^ 

~ 

disintegration as illustrated in - 
figure 3.29, $T 

o 

Zoge_2^ This zone is a strip 

on the roof of width equal to one- o 
quarter of the length of the struc- 
ture measured from the edge of Zone 
1 toward the center of the roof. o 

Figure 3,32 is a plot of the 
ratio P roof /P s for any point a 
distance L f from the front edge 



Figure 3*31 is a plot of 

-.-n^vpf-eJ the ratio of the local 
roox s 

overpressure at any point on the 
roof in terras of time t to the 
overpressure in the incident blast 

wave where t, L f /U is the time- 

d ' o 

displacement factor > equal to the 
time required for the shock front to 
travel a distance L f , the distance 
from the front edge of the roof to 
the point under consideration. 

Figure 3*31 indicates that the 
local overpressure versus time curve 



Bias! Wov@ 



C.L. 



P'may not 
b !ts8 than 0.5. 




Timt.t 

Figure 3.32. Local roof or side wall over- 
pressure ratio vs time 9 Zone 2~closed 
rectangular structure 
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Here again, 
t. 



incorporate the time- 



L f /U . 
1 o 



P . and P 
roof s 

t is given by equation (3-ll). 



of the roof, 
displacement , 

Zone 3* This zone includes all points on the roof not included in 



roof' 



for all 



H-H 



Note! 



Zones 1 and 2. Figure 3*33 is a plot of the ratio 
points in this zone. Zone 3 is 
essentially a region of two- 
dimensional shock phenomena for 
which considerable data are avail- 
able [ 4, 5, 9, 10, 11, 12, and 13]. 

The effects of the vortex 
travel across the roof of the 
structure are confined to Zones 
2 and 3. The difference in the 
overpressure existing in these 
zones is due to the assumed 
severity of the vortex effects* 

As the shock front passes 
over the point being considered, 
the local overpressure is raised 
to the overpressure in the inci- 
dent blast wave. This equality between local and incident blast wave over- 
pressures is maintained until the vortex developed at the front edge of the 
roof detaches itself and moves toward the rear edge, causing a decrease 
from the overpressure in the incident blast wave at the point being con- 
sidered. The local overpressure reaches its minimum value at the time that 
the vortex is over the point in question. The time required for the vortex 

to travel the distance L f is 

(3.H) 




Time, t 

Figure 3.33* Local roof or side wall over- 
pressure ratio vs time. Zone 3 closed 

recKmgiilar structure 



t m L'/v 



where v is the vortex travel velocity, obtained from a reduction of the 
data in reference [^-] and given by the relation 



fo. 



Olj-2 



0.108 - 



U 



The value of the ratio of P 



rooi 



p 



Q 
at time 



m 



is 



(3-12) 



(3.13) 
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For points located in Zone 2 the minimum value of equation (313) is 
0.5; for points located in Zone 3 the minimum value is zero. 

As soon as the vortex has passed the point in question, the local 
overpressure starts to build up until at time t t m + 15h ! /U Q (where h 1 
is the clearing height of the structure) it is once again equal to the 
overpressure in the incident Shockwave. 

e. Average Roof Overpressure P -. Having established the varia- 
tion of local overpressure, the average overpressure on the roof at any 
time could be obtained by the summation of the local overpressure curves 
over the entire roof at a given time t. If there were no lateral varia- 
tion of local overpressures, such a method could be used to develop a 
general procedure applicable to all structures. However, the presence of 
this lateral variation complicates any procedure to such a degree that oiily 
the limiting case for Zone 3 loading is considered. The ratio of the aver- 
age roof overpressure P^Q^ " bo 
the overpressure of the incident 
blast wave existing at the 

center of the roof P is 

s 

plotted in figure 3.3^ for times 

in excess of time 

where 



NOTE: P"0.9+0.l(l.O--~2); except P" may not exceed 1.0 

P 1 * 2jO-(r~~ + 1 XT") or 0.5+ 0.125 (2- ~z) 
whichever i* smaller; except P 1 may not be lest 
than zero. 

For 0t jL/U , P r oof n d ^*id* ***** 
linearly from zero at t*0, to the value 
indicated at t > L/U . 




t - L/U o 
P incorporates the 

L/2U o . 



time -displacement t. 



The average overpressure on the 
roof varies linearly from zero 
at time t to the value of 
P roof given in figure 3-3^ at 



Time, t 

Figure 3.34. Average roof or side wall 

overpressure ratio vs time, Zone 3 

closed rectangular structure 



time 
that 



L/U . The assumption 



roof 



can be expressed as 



a proportion of P is valid 

s 



only if the time required for 
the shock front to travel the length of the building is small compared to 
the duration of the positive phase. For values of (L/U Q )/t o less than 0.1 
this assumption is reasonable. ^~ 

A farther restriction which must be imposed is due to the importance 
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of the lateral variation of overpressures on the roof. If the width of the 
structure normal to the direction of travel of the shock wave is greater 
than twice its length, the average roof overpressures as determined by 
figure 3*3^ are satisfactory. For a structure whose width is less than its 
length, the average roof overpressure for times greater than time t L/U 
is more correctly given by the relation 

P - * P (3-1^) 

roof s w / 

where P incorporates the time -displacement t , L/2U . The average 
roof overpressure varies linearly from time t * to the value calculated 
by equation (31*0 sit time t L/U . For structures which are approxi- 
mately square in roof plan, neither method is the more correct and an aver- 
age of the two should be used. If a more accurate value of the average 
overpressure on a portion of the roof, such as a roof slab, is desired, 
this may be obtained by computing the local overpressure at two or three 
points on the slab and obtaining the weighted average as a function of 
their common time t. 

f. Local and Average Side Wall Overpressure P ., and P ., On 
_ _ _ f _ side _ side. 

the basis of test data it has been observed that the sides of a structure 
are loaded in the same manner as is the roof, and that figures 3-31 to 3*3^- 
determined for the roof of a structure apply equally well to the sides. 

The side walls are divided into three zones,* as shown in figure 3*30* 
for the determination of local overpressures P .. . Figures 3-31 to 3*33 
present these local overpressures in the form P . . /P , where P and 
P incorporate the time-displacement t L f /lJ . L ! is the distance from 
the front edge of the side wall to the point under consideration. 

The average overpressure on the side is given as the ratio P /P 

in figure 3-3^ where P is the overpressure in the incident blast wave 

s 

with a time-displacement t, L/2U . This figure is applicable only if 
the height h of the side wall is greater than half the length L. If 
this restriction is not satisfied the average overpressure on the side is 
more correctly given by the relation 



3-09 PROCEDURE FOR COMPUTATION OF LOADS ON CLOSED RECTMGUIAR STRUCTURES. 
The following paragraphs present a series of step-by-step general procedures 
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for the computation of the "blast wave overpressures on the various por- 
tions of a closed rectangular structure. 

Step 1. Given data: Size of weapon j distance of proposed structure 
from ground zero; exterior dimensions of structures. 

Step 2. Determine P for the given data from figure 3.12a or 
3*12b using the scaling relations of paragraph 3-06 as necessary. 

Step 3* Determine t . knowing P , by applying the scaling rela- 

o so 

tion of equation (3*3) to figure 3.10a or 3*10b. 

Note that all the quantities listed in Step 1 are not necessary to 

determine the loading on any structure if P and t plus the exterior 

so o 

dimensions of the structure are given. 

a. Average Front Wall Overpressure vs Time. Step 1. Determine P 

111 w * -a- 3 

for a sequence of times t.., tp, t,-, etc., from table 3.1 or figure 3-22 

for t, m 0. 
a 

Step 2. Determine P ., knowing P , from figure 3*20. 

re i j_ so 

Step 3. Determine the clearing height of the front wall h 1 . This- 
is equal to the full height of the front wall or half its width, whichever 
is the smaller. 

Step 4. Determine c -, knowing P , from figure 3.21. 

re i j. so 

Step 5- Determine t from the relation t 3/h ! c -. 

c c rei JL 

Step 6. Determine q Q , knowing P , from figure 3,23. 
Step 7- Determine q for a sequence of times t., tp, t_, etc., 
from table 3.2 or figure 3.2^. 

Step 8. Determine (P + 0.85q) for the sequence of times used in 

s 

Steps 1 and J above. The following headings systematize the calculations 
of Steps 1 to 8. 

time, t t/t o q/^ q P^ Pg Pg + .85q 

(1) (2) (3) (*) (5) (6) (7) 

Column (l) is a set of times ranging from t to t t suffi- 

o 

cient in number to allow drawing of the curve of P + 0.85q. 

s 

Column (3) is values obtained from table 3.2 or figure 3.24 for the 
values of t/t Q ift column (2). 

Column (5) is values obtained from table 3.1 or figure 3.22 for the 
values of t/t in column (2). 



I 
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Step 9. Plot the curve of P front vs time as illustrated in figure 
3.25. 

t>. Average Back Wall Overpressure vs Tine. Step 1. Determine TJ , 
knowing P , from figure 3.9. 

SO 

Step 2, Determine t,. from the relation t^, L/U . 

Q. d ' o 

Step 3. Determine t. from the relation t, * kh'/c . where c 

D D ' O' O 

1115 fps and h f is either the full height of the back wall or half its 
length, whichever is the smaller. 

Step k* Determine P , at time t - t, * t, from table 3*1 or 
figure 3-22. 

Step 5* Determine (P- back ) mx at time "t * t + t from the rela- 

tion ^ s I P sb 



The average overpressure on the rear face varies linearly from zero 

at time t . t, to (P, _ ) at time t t^ + t . 
d % back'max _ d TD 

Step 6* Determine the values of P, , /P for a sequence of times 



in excess of t t , + t from the relation 



, 



p (p ) 

back _ back max 

P " P , 



1 - 



(p ) 

x 



back max 



t - 



t - 

o 



Step 7. Determine the values of P for the sequence of times in 

s 

Step 6^ and hence the values of P, , . The following tabular headings 
systematize the work in Steps 6 and 7- 

time, t t - t , (t - tj/t P /P P EU ,/P PU i 
' d x d y/ o s 7 so s back 7 s back 

(D (2) (3) (10 (5) (6) (?) 

Column (l) is a sequence of times ranging from t = t, + t, to 

t t + t, . 
o d 

Column (k) is a set of values obtained from table 3.1 or figure 3.22 
for the values listed in column (3)- 

Column (6) is the values calculated from the relation in Step 6 for 
the times given in column (l). 

Step 8. Plot the curve of P, , vs time as illustrated in figure 
. 

c. Average Net Horizontal Overpressure vs Time, Step 1. Knowing 
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P f ont and P ba k for a se< l uence o:f ^mes *-,* "tp' ^V etc., determine * 

5 net = *Wt ' *ba<* fOr each ^ ue of *' 

Step 2. Plot the curve of P , vs time as illustrated in figure 

net 

3.28. ; 

d. Local Roof Overpressure vs Time, The step-by-step procedure 
listed belov is for a point located within Zone 3 (see figure 3-30) of the 
roof. The procedure for Zone 2 is similar except that figure 3.32 applies 

instead of figure 3*33, and the value of the ratio P /P at time 

roox s 

t * t (Step k) may not be less than 0.5. The procedure for Zone 1 in- 
volves a simple computation of P at various times using table 3.1 or 

s 

figure. 3. 22 as illustrated in figure 3-31. 

Step 1* Determine t, - L'/U - 

do 

Step 2. Determine the vortex velocity v from the relation 



( 0.01*2 + 0.108 ~ ]U 
\ L / o 



Step 3. Determine t = L r /v. 

m ' 

Step k. Determine the value of the ratio P _/P at time t * t 

/- t \ roof s m 

from the relation P roof /P s = MP BO A^-T)l^ - !) + 1-0, except that the 

mini muni permissible value of this ratio is zero. The ratio P /P de- 
roof' s 

creases linearly from 1.0 at time t (t, + t )/2 to P Jp at time 

d m ' roof' s 

t and increases linearly to 1.0 at time t t + 15h f /U i for all other 
m ni o 

times the value of the ratio is 1.0. 

Step 5. Determine the value of the ratio P /P for a sequence 

rooz s 

of times between time t * t , and time t t, + t . 

d do 

Step 6. Determine P g for the sequence of times in Step 5. 

Step 7. Determine P roof from Steps 5 and 6. The following tabular 

headings systematize the vork in evaluating P 

roof 

time, t t - t (t - t,)t P /P P P /P p 

d d x o s- so s roof 7 s roof 

(D (2) (3) (k) (5) (6) ( 7 ) 

Column (l) is a sequence of times ranging from t t to 

d 

t t . + t . 
a o 

Column (k) contains the values obtained from table 3.1 or figure 
3-22 for the values in column (3). 

Column (6) contains the values obtained in Step If above. 

46 
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Step 8. Plot the curve of P rQof V s time from the values in column 
(7) as illustrated in figure 3.33. 

e Average Roof Overpressure vs Time. The step-by-step procedure 
listed below is for a structure for which figure 3.3^ is applicable, as 
explained in paragraph 3~08e. 

Step 1. Determine t, L/2U . 

Step 2. Determine the value of P /P at time t t, + L/2U n 

roo s^ 

from the relation P roof /P g - 0.9 + 0.1 (l.O - l <1.0 



Step 3. Determine the value of P _/P at time t 5L/U from 

roof s ' o 

the equations P/P . 2.0 - fe+ l feY * or P/P . 0. 5 + 



O 

I P so \ 
0.125 (2 - ilj.,7/ whichever gives the smallest value, but not less than 



zero. 



l 



Step k. At time t 5L/U + 15h f /U the value of P Jp is 
_ o x 7 o roor s 

one. P oof / p s varies linearly between the values computed in Steps 2, 
3, and U, and remains equal to one for times greater than t 5L/U + 
15h'/U Q . 

Step 5- Determine the values of P /P for a sequence of times 

rooi s 

greater than t L/U by interpolating between the values computed above, 

Step 6. Determine P for the sequence of times used in Step 5- 

s 

Step 7- Determine P from Steps 5 and 6. The following tabular 
headings systematize the work in evaluating P f . 

time, t t - t d (t - t d )/t o P s /P sQ P s P rOQf /P s P roof 
(1) (2) (3) (V (5) (6) (7) 



Column (l) is a sequence of times ranging from time t L/U to 

time t = t + t , . 
o d 

Column (k) contains the values obtained from table 3.1 or figure 3.22 
for the values in column (3)- 

Column (6) contains the values obtained from Step 5 for the times in 
column (l) . 
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Step 8. Plot the curve of 



3-09* 



vs time as illustrated in figure 



3.34. Note that P 



varies linearly from zero at time t to the 



L/U Q . 



roof 
computed value at time t 

f . Local Side Wall Overpressure vs Time. The local side wall over- 
pressure- for a point in Zone 3 of the side is computed by use of Steps '1 to 
8 of paragraph 3 -09d. The procedure for Zone 2 is similar except that 
figure 3.32 applies instead of figure 3-33> and the value of the ratio 

P /P. at time t * t (Step if) must not be less than 0.5* The proce- 
roox s m 

dure for Zone 1 involves a simple computation of P at various times 

s 

using table 3-1 or figure 3.22 as illustrated in figure 3-31. 

g. Average Side Wall Overpressure vs Time. The average side wall 
overpressure for a structure for which figure 3-3^ is applicable is com- 
puted by use of Steps 1 to 8 of paragraph 3~09e. 

3-10 NUMERICAL EXAMPLE OF COMPUTATIONS OF LOADS ON CLOSED RECOT^GUIAR 



STRUCTURE. 



Blast 
Wave 



27ft. 



^-Front Wall 



Roof-" 



Back 



15ft. 



-54ft.- 



Ground 



Data: 18 KT typical airburst. 

3000 ft distance from ground 
zero. 

Exterior dimensions of prob- 
lem structure (figure 3*35). 



Note: Length of Bldg. Parallel to Shock 
Front is 200 ft. 

Figure 3.35. Exterior proportions closed 
rectangular structure 



From figure 3.12 P 



so 



Prom figure 3.10a t 
(1 ET for 10 psl). 01 

For an 18-KT weapon, t 
(0.262) * 0.685 sec. 



10 psi. 
* 0.262 



18' 



1/3 



a. Average Front Wall Overpressure vs Time. Step 1. For P 
10 psi from figure 3-20, P refl 25*3 psi. 

Step 2. From figure 3-35, h f * 15 ft. 

Step 3. For P * 10 psi from figure 3.21. c _. = 1290 fps. 
so refl ^ 

Step k. t c 3h'/c refl = 3(15)/1290 - 0.03^9 sec. 

Step 5. For P gQ 10 psi, from figure 3.23, <J Q 2.23. 
Steps 6 through 8. In table folloving. 



so 



I 



3-10b 





0.0685 
0.137 
0.2055 
0.27^ 
0.3^25 



EM 



0.1*795 



0.6165 
0.685 





0*10 
0.20 
0.30 
O.UO 
0.50 
0.60 
0.70 
0.80 
0.90 
1.00 



2.23 



0.89 
0.55 
0.33 
0.19 
0.11 

0.06 

0.03 
0.01 
0.00 



P s/ P so 

1.000 
0.81U 

0.655 
0.519 

0.1KJ2 
0.303 
0.220 



O.OQO 

o.oin 

0.000 
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p + o.85a 

s 
12.12 



71 
33 
21 
30 
1.55 
0.93 



0.00 



0.00 



Sample calculation of typical 

m 0.2714.. 

Prom table 3.2 for t/t Q 
O.UO, q/q. - 0.1U8 
q, = 0.lW(2.23) - 0.33 psi 

From table 3*1 for t/t Q 
. O.to, P/P - 0.1*02 



entry in table: For t/t Q 0.1*0(0.685) 



SSQ 
P = 0.^)2(10) 1^.02 psl 

P s = 0.85q 



02 + 0.85(0.33) 

- 14-- 33 psl 

Step 9* Plot the curve as 
shown in figure 3.36(a)- 

b. Average Back Wall Over- 
pressure vs Time. Step 1. For 
P = 10 psi, from figure 3.9, U Q 



t - L/U 



1^03 fps. 

Step 2.. 
0.0385 sec. 

Step 3- 
U(l5)/ni5 * 0.538 sec. 

Step 

0.0538/0.685 = 0.0785, from 
0.852, 



For (t - t d )/t Q 



table 3.1, P g /P so 




O.I 03 0.5 0.7 

Time (tec) 

(a) Average Overpressure 
vs Time Curve 



O.I 0.3 0.5 0.7 

Time (sec) 

(b) Aver age Net Horizontal 
Overpressure vs Time 
Curve 
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x 



r sb 



10(0.852) 8.52 psi. 



O.I 0.3 0.5 0.7 

Time (sec) Time (sec) 

(c) Local Roof Overpressure (d) Average Roof 

vs Time Curve Overpressure vs Time 

Curve 



Figure 3.36. Overpressure vs time curves-- 
closed rectangular structure 
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Step 5- ( 



back'max 



Steps 6 and 7- 



3-10c 

P gb [l + (1 - P)e~ P ]/2 where P 0.5P so /l4.7. 
p = 0.5(lO)/l4.7 0.34 

I = 8.52 [l + 0.47] /2 - 6.26 
max 

In table following. 



t 


t - t d 


V* " % dJ' 


V r so 


*s 


back' s 


back 


0.0923 
0.1070 

0.1755 
0.2440 
0.3125 


0.0538 
0.0685 

0.13T 
0.2055 
0.274 


0.079 
0.1 
0.2 

0.3 
0.4 


0.850 
0.814 
0.655 
0.519 
0.402 


8.50 
8.14 
6.55 
5.19 
4.02 


0-735 
0.735 
0.740 
0.750 
0.767 


6.26 
5-98 
4.85 
3.89 
3.08 


0.3810 
0.4495 


0.3425 
o.4n 


0.5 
0.6 


0.303 
0.220 


3.03 
2.20 


0.790 
0.820 


2.39 
1.80 


0.5180 
0.5865 
0.6550 


0.4795 
0.548 
0.6165 


0.7 
0.8 
0.9 


O.X49 

0.090 
0.041 


1.49 
0.90 
0.41 


0.855 
0.897 
0.947 


1.27 
0.81 
0.39 


0.7235 


0.685 


1.0 








1.000 






Sample calculation of typical entry in table above: 

For (t - t d )/t Q = 0.5, t = 0.5(0.685) + 0.0385 = 0.3810 
Prom table 3-1, for (t - tj/t = 0.5, ?./?. = 0.303 



P S = 0.303(10) = 3-03 psi 
For (t - t d )/t o =0.5 



"bac 




s' so 



' (P back W P sb 




^1 [0.3810 - (0.0385 + 0.0538)' 
''* J L 0.685 - 0.053S 



)1' 
J 



= 6.26/8.52 +[l - 6.26/8 

= 6.26/8.52 + [l - 6.26/8.52] [0. 2887/0. 6312] 2 

= 0.790 

p back - 3-03(0.79) = 2.39 Psi 

Step 8. Plot the curve as shown in figure 336(a). 
c. Average Net Horizontal Overpressure vs Time. Step 1. Determine 

values of P. , and P, . at times t = 0.00, 0.05, 0.15, 0.20, ..., 
front oacK _ 

0.65, 0.70, from figure 3-36(a) and determine P. _^ = ?____,. - ?,_._,_ foir 
each time. 



P = p - p 
net front 



Step 2. Plot curve of P vs time as shown in figure 3.36(t>). 

net} 

d. Local Roof Overpressure vs Time at Point x of Figure 3.35. 



Step 1. t, = 



= 27/llK)3 = 0.0192 sec. 
50 
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Step 2. v = (o.Ote + 0.108 ~)U Q = [0.0^2 + 0.108 (0.5)] 1^03 

= 13^.7 fps. 

Step 3* t = L'/v = 27/13^. 7 = 0.200 sec. 
Step h. A? time t t^ P rQof /P s - MP flO /3A.7) (^ - l) + 1 



rQofs flO 

p roof /p s = K 10/1^.7 )(27M - 1) + i = -0.36 < o. 

Therefore. P JP =0 at time t = 0.200 
7 roof s 

P Jp = 1.00 at time t = t + IJh'/lJ 
roof' s m _ o 

Steps 5 through 8. In tahle following. (Note that in the computa- 
tions in table all times in columns (l) and (2) are calculated from assumed 
values in column (3) except at times t = t m = 0.200, t = T> (* d + t m ) = ^ 
(0.0192 + 0,200) = 0.1096, and t = 0.3602.) 



t 


*'*! 


(t - t d )/t 


P s/ P o 


P s 


P /P 
roof 7 s 


P 
roof 


0.0192 








1.000 


10.00 


1.0 


10.00 


0.08TT 


0.0685 


0.1 


0.814 


8.14 


1.0 


8.14 


0.1096 


0.0904 


0.132 


0.761 


7.61 


1.0 


7.61 


0.1562 


0.137 


0.2 


0.655 


6.55 


0.483 


3.16 


0.200 


- 


- 


- 


- 








0.22l)-7 


0.2055 


0.3 


0.519 


5-19 


0.154 


0.80 


0.2932 


0.274 


0.4 


0.402 


4.02 


0.583 


2.34 


0.3602 


0.3410 


0.498 


0.305 


3.05 


1.0 


3-05 


0.3617 


0.3425 


0.5 


0.303 


3.03 


1.0 


3-03 


OA302 


0.411 


0.6 


0.220 


2.20 


1.0 


2.20 


0.4987 


0.4795 


0.7 


O.l4"9 


1.49 


1.0 


1.49 


0.5672 


0.548 


0.8 


0.090 


0.90 


1.0 


0.90 


0.6357 


0.6165 


0.9 


0.041 


0.41 


1.0 


0.41 


0.7042 


0.685 


1.0 








1.0 






Sample calculation of typical entry in table above: 

For (t - t d )A = 0.4, t = 0.4(0.685) + 0.0192 = 0.2932 sec 
From table 3.1, for (t - t,)/t = 0.4, P /P cn = 0.402 

OL O S SO 

P = OA02(10.0) = k.02 psi 
s 

Since P _/P varies linearly from at t = 0.200 to 1.0 at 
roof' s 

t = 0.3602 
P roof/ P s = (' 2932 " 0-200)/(0.360 - 0.200) = 0.583 



P roof = -583(4.02) = 2.34 psi 

Step 7. Plot the curve as shown in figure 336(c). 
e. Average Roof Overpressure vs Time. Step 1. t^ L/2U Q 
54/2(1403) 0.0194 = 0.0194 sec. 

Step 2. At time t - L/2U Q + t d -- 0.0388 sec 

51 
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P /P . 0.9 + 0.1(1.0 - P sn A^7) 2 - 0.9 + 0.1(1.0 - 10/14.7) 2 - 0.910 
roor s \ 

(i.e., this value may not exceed one;. 

Step 3. At time t 5L/U Q - 0.1928 sec. Compute 2.0 - 
L + !j|SL\ [|1J 1/3 2.0 - (l + ^\ (15M) 1/3 - 0.902 and 0.5 + 

V / \ * f* 



0.125 2 - 



0.5 + 0.125 2 - 



\ ' 

Stoce 0.717 < O'.90e, Ppoof/Pg = 

Step 4. At tine t - 5L/U Q + 




0.717 



15h'/U o - 



,1928 + 0.1604 - 0.3532 sec 



*W/ P s ' X ' 

Steps 5 through 7. 



In table following. 



t 


t - t d 


(t - t d )/t o 


p /p 
s' so 


P s 


p /p 
roof 7 s 


P roof 


0.0388 

0.0879 
0.1176 


0.0194 
0.0685 
0.137 


0.0283 
0.1 
0.2 


0.944 
0.814 
0.655 


9.44 
8.14 
6.55 


0.910 
0.848 
0.762 


8.58 
6.90 
4.99 


0.1928 


0.1734 


0.252 


0.581 


5.81 


0.717 


4.17 

i .^ _ 


0.2249 
0.2934 


0.2055 
0.274 


0.3 
0.4 


0,519 
0.402 


5-19 
4.02 


0.774 
0.894 


4.01 
3.60 


0.3619 


0.3425 


0.5 


0.303 


3.03 


1.0 


3.03 


0.4304 


o.4n 


0.6 


0.220 


2.20 


1.0 


2.20 


0*1^989 


0.4795 


0.7 


0.149 


1.49 


1.0 


1.49 


0.5674 


0.548 


0.8 


0.090 


0.90 


1.0 


0.90 


0.6359 


0.6165 


0.9 


0.041 


0.41 


1.0 


0.41 


0.7044 


0.685 


1.0 








1.0 






Sajrfple calculation of typical entry in table above: For (t - t d )/t Q 
= 0.1, t 0.1(0.685) + 0.019^ - 0.0879 sec 

From table 3-1 for (t - t d )/t Q 0.1, P G /P 



s' so 



.1U psi 



Since P /P varies linearly from 0.910 at t 0.0388 sec to 
roox s 

0.717 at t = 0.1928 sec 

P /P - Q10 (Q-OQ79 - 0-0388 \ , 
P roof /P s ~ ' 91 " \ 0.1928 - 0.0388,/ ( ' 

= 0.910 - O.o49l(0.193)/0.l54o 
= 0.848 
P~roof = 0.848(8.14) =6.90 psi 

Step 8. Plot the curve as shown in figure 3.36(d). 
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LOADING ON KECTANGUIAR STRUCTURES WITH OPENINGS 



3-11 EFFECTS OF OPENINGS ON LOADING, For the determination of the loads 
on rectangular structures with openings it is necessary to consider all 
structures as falling into one of two categories. 

Category (l). Structures having open interiors free of walls and 
other obstructions so as to allow relatively unhindered propagation of the 
"blast wave through the interior of the structure, once it has entered 
through windows or otherwise* It is assumed that the percentage of "back 
face openings is approximately equal to the percentage of front face 
openings . 

Category (2). Structures having interior walls and other construc- 
tion which would prevent free propagation of the blast wave through the 
interior of the structure. 

Those structures falling into Category (2), above, are analyzed in 
a ifianner similar to that for the case of closed rectangular structures when 
determining the net trans lational forces acting on the structure. Net 
loads on exterior and interior portions of these structures may "be computed 
as are the loads for structures falling into Category (l), with suitable 
interpretation of the various structural dimensions involved. The deter- 
mination of the loads on structures with relatively open interiors, Cate- 
gory (l), is described below for approximately equal percentages of open- 
ings in the front and. "back walls . 

The sequence of phenomena following the impingement of a shock wave 
on a rectangular structure with openings can be briefly summarized in the 
following paragraphs. 

When the shock front strikes and reflects from the front wall, the 
overpressure on that wall rises *to the reflected overpressure. Rarefaction 
waves Immediately move in from the outside edges of the wall, clearing the 
reflected overpressures. 

Windows and doors in the front wall will probably break under the re- 
flected overpressure, though not instantaneously. However, they generally 
will break before the clearing is complete; hence, part of the clearing 
will occur through the window openings, allowing high pressure air to flow 
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into the Interior. This sudden release of high pressure will cause shocks 
to form in the interior. These inside shocks, formed at each opening, will 
spread downstream from each opening and will tend to combine into a single 
shock, front* This interior shock is, initially, weaker than the incident 
shock* 

Meantime the incident shock front has moved past the front wall and 
is sweeping across the exterior roof and sides. The shock then moves 
across the exterior of the back wall of the structure from the sides and 
from above, building up the overpressure on that surface. 

For those structures which have no interior partitions to obstruct 
the passage of the shock through the interior, the shocks entering at the 
front openings will move through the structure raising the overpressures on 
the inside surfaces as they are covered. 



a. Average Exterior Front Wall Overpressure P 



front* 



The time re- 



quired for reflected overpressures to clear on the front of a solid wall 
struck by a shock wave is expressed in multiples of the time necessary for 

a rarefaction wave to sweep the wall once, h'/c .. When walls with 

rex JL 

numerous openings are considered, however, clearing can take place around 
the edges of the openings. This necessitates a modification of the method 
of determining clearing times because the rarefaction wave travels a much 
shorter distance to cover the wall. 

The distance hi is introduced as the weighted average distance the 
rarefaction wave must travel to cover the wall once, assuming Jbnmediate 

u * j access of the incident shock to the interior 

of the structure. As shown in figure 3.37 
the front face is divided into rectangular 
areas, determined by the location and dimen- 
sions of openings considering the directions 
along which the axea can be cleared in the 
shortest possible time. These areas are 
labeled in the following manner: 

Step 1. 
from two opposite sides. 



Opening 



Opening 



^ 



Figure 3.37. Subdivision of a 

typical wall with openings to 

determine weighted average 

clearing height 



Label 1 , all areas cleared 



Step 2. label 



all areas cleared from two adjacent sides. 
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Step 3 Label 2, all areas cleared from one side. 
Step 4. Label 3* all remaining areas. 
The weighted average clearing height is then: 

8 h A 

* 1 '' (3.16) 

where A f net area of front face (wh less openings) 

A .a area of each portion of the subdivided front face, except 
openings 

h = for areas 1 , the average distance "between the sides from 
which clearing occurs 

* for areas 1, and 3, the average height or width, whichever is 
smaller 

= for areas 2, the average distance between the side from which 
clearing occurs and the side opposite 

5 a the area clearing factor. 

1/2 for areas 1 
1 a 

* 1 for areas L 
1 for areas 2 

* 1 for areas 3 6 h A 

Z the notation to represent the summation of the ' n n n 
term for all areas f 

h f = h or tf/2, whichever is the smaller 

The time required to clear the front face of a structure with open- 
ings of reflection effects t ! can now be evaluated from the relation 

*i - 3h f/ C refl ^3.17) 

Windows in a structure .will not be immediately blown out; GEREEMHOUSE 
data (see page 72, reference 15) indicates that the window-breaking time is 
of the order of a millisecond or less. This time, however, is too short to 
be of any practical importance in computing the loading on the structure. 
Therefore, it is assumed that the windows fail as soon as they are struck 
by the shock front and do not affect the loading computation. 

The time c\irve of the average outside front wall overpressure on the 
net wall area can be constructed by the methods outlined in paragraph 3-08a 
and illustrated in figure 3*25 using the clearing time as defined by 
equation (3.17). 
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b. Average Interior Front Wall Overpressure P* .- + ^he average 

interior overpressure on the inside of the front vail P initially 
builds up in a similar manner as the average overpressures on the exterior 
back vail of a closed structure, except the vortices are located all around 
the inside edges of all openings* The build-up is linear from zero at time 
t to the peak overpressure P . (see paragraph 3-He and figures 
3-^1 and 3A2) of the interior shock at time t ^h' /c > vhere h' is 

the veighted average build-up 
height of the interior front 
vail computed from equation 
(3.16) but vith all quantities 
interpreted as applying to the 
interior surface of the front 
vail. The overpressure is 
maintained aiid builds up again 
vhen the interior shock is re- 
flected back from the rear 




Time, t 



Figure 3.38. Average interior front wall overpressure 
vs time small rectangular structure (L./C io ) < 0.1t Q 



vail to a peak value vhich is 
dependent on the size of the 
structure. Equilibrium is 



then attained vhen the overpressure vithin the interior reaches P ; P 

s s 

being the blast vave over- 
pressure existing on the 
outside of the structure. 
Figure 3.38 illustrates the 
overpressure vs time curve 
for the inside front face 
of those structures in vhich 



10. 

r 

3 

w 

M 



io 



< 0.1 t 



and 



figure 3.39 for those in 



vhich 



> 0.1 t 




Figure 3.39. Average interior front wall overpressure vs 
time large rectangular structure (L./li. o )> 0.1 t Q 



c ' Average Net Front 
Wall Overpressure P- . . The average net overpressure acting on the 



I 



3-Hd 



front wall P, 
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is deter- 



"f-net 
mined by subtracting the interior 

front wall overpressure from the 
exterior front wall overpressure 
on a common time basis according 

to the equation 

p . p - p 
f-net front i-fr 

where P , , P ,, and 

PJ t are given in terms of 

t. figure 3-^0 illustrates 



f-net 



= P 



front 



_ p 



i- front 




Average Exterior Front Wo 1 1 Overpressure, P fronf 
Interior Front Wall Overpressure, P { _ front 

AverogeNet Front Wall Overpressure,?* 

f ~ ntt 



Time, t 



graphically the variation of the 
average net front wall overpres- 
sure with time. 

d. Average Exterior Back Wall Overpressure P, 



Figure 3.40. Average net front wall overpressure 
vs timerectangular structures with openings 



A time lag is 



present between build-up of pressure on the back wall due to the time re- 
quired for the blast wave to travel the length of the structure. Until the 
windows or window frames in the back wall fail, the back wall loads behave 
in the same manner as those on the back wall of a closed structure. After 
the windows break the back wall loading is affected by the vortices formed 
at all of the edges of the rear openings due to the interior blast passing 
through the openings at the same time that the exterior wave sweeps the 
rear face. If it is assumed that the windows blow out instantaneously over 
the entire structure, i.e., the window-breaking times are zero, the build- 
up time for pressures over the rear face is 



where 



is the weighted average build-up time of the rear face, defined 

5 h A 



(3.20) 



Ap is the net area of the back wall (gross area less openings) 
8 , h , A are as defined in paragraph 3~11& 

h/ is computed in the same manner as hi. (the weighted average 
clearing height of the front face) 

The average overpressure vs time curve for the outside rear face can 
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be computed by the methods outlined in paragraph 3~08b and illustrated in 
figure 3-27 using the build-up times as given in equation (3*19)- 

e. Average Interior Back Wall Overpressure P-.-u The presence 
of openings in the front vail of a structure permits a portion of the inci- 
dent shock to enter through these openings, producing an interior shock 
wave veaker than the incident, which becomes approximately plane as it 
travels the length of the interior. That portion of the interior shock not 
passing through the back vail openings is reflected by the back vail, ac- 
cording to the lavs of normal reflection and returns tovard the front. 
Quick clearing occurs on the inside and vortices form at edges of openings 
on the outside as the blast vave passes out through openings. 

In the prediction of inside pressures the structure is assumed to be 
smooth* Structural features vhich do not follow this idealization are 
neglected as regards their influence on the determination of interior 
loads. These features include interior columns along side vails, pipe 
sections vhich cross the interior, etc. The peak overpressure of the in- 
terior shock wave P is dependent upon the overpressure of the incident 
shock vave and the percentage of vindov openings. Figure 3.4l is a plot of 
P goi vs P SQ for different values of the ratio of the area of openings in 

the front face to the gross area of the front face A JA. _,. Figure 3.^2 

of' gf 

is a plot of P gQi vs A of /A for various peak incident overpressures. 
These curves are based upon data reduced by the Armour Research Foundation 
[l7] from shock tube tests performed at the University of Michigan [l6] and 
at Princeton University [l8] . 

The inside back vail remains unloaded until time t * L /U , vhere 

d i' 10* 

1^ is the distance from the outside of the front wall to the inside rear 
face and U iQ is the velocity of the inside shock front vhich is plotted 
in figure 3.9. The overpressure on this vail is then instantaneously 
raised to the reflected overpressure value P given in figure 3.20. 

L""".L 6X JL 

The presence of openings in the back vail prevents the maintenance of 
this reflected overpressure and, hence, there is a linear decay to a value 
approaching the drag pressure in the inside blast vave. The time required 
for this decay is a function of the Dumber and size of openings in the back 
vail and of the size of the structure. For relatively small structures 
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there is an additional build-up following an initial decay due to the re- 
turn of the shock wave reflected from the rear wall toward the front wall, 
while for longer structures this second build-up is absent. It is assumed 
that the interior blast wave eventually reaches a condition of equilibrium 
wherein the interior and 
the incident blast waves 
have equal overpressures p 
and the back wall inte- g 

JO 

P 1 



rior overpressure equals 
the drag pressure in the 
incident blast wave, with 
the time displacement 
factor t d * L i/ U i - 

To account for the 
effect of size of the 
structure on interior 
overpressures, two curves 
of the interior back wall 
overpressures are given. 
Figure 3 A3 is applicable 
to structures in which the 



|OL 
T 

P 80 +0.85q 



O 




C 



Time, t 




Figure 3.43. Average interior back wall overpressure vs 
time small rectangular structure (L i /V io )< O.lt o 



time required for the inside 
reflected shock wave to 
travel the length of the in- 
terior is less than one- 
tenth the duration of the 
positive phase of the out- 
side shock wave; i.e., 

L i/ U io < - lt ' Fi sure 3A1+ 
is applicable to those struc- 
tures for which /^ > 



C 



O.lt 



as a func- 



Figure 3.44. Average interior back wall overpressure vs 
time large rectangular structure (L./U . ) > O.lt 



i-back 

tion of time is as illus- 
trated in these figures. 
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NOTE: P b -ntt = Pi- back - P bock 



I 



O 




3-Uf 

f . Average Net 
Back Wall Overpres- 



sure 



b-net * 



The 



Average Interior Back Wall Overpressure, Pi-bock 



average net overpres- 
sure acting on the 



Average Exterior Back Wall Overpressure, P bock 
X/Average Net Back Wail Overpressure, P b-net 



back wall 



is 



Time, t 

Figure 3.45. Average net back wall overpressure vs time- 
rectangular structure with openings 



b-net 



P i-back " P back 



b-net 
determined by sub- 
tracting the exterior 
back vail overpres- 
sure from the inte- 
rior back wall over- 
pressure on a common 
time basis according 
to the equation 

(3.21) 



where P b _ net , ^.^0^ and S^ck are STen in tems or t * F1 S ure 3- 1 
illustrates graphically the variation of the average net back wall over- 
pressure with time. 



g. Average Exterior Roof Overpressure P 



roof* 



In this paragraph it 



is assumed that the roof is plane, horizontal, and vithout roof surface 
openings. Therefore, the average overpressure vs time curves on the ex- 
terior are computed by the methods outlined in paragraph 3~08e and illus- 
trated in figure 3-3^ The quantity h 1 used in these figures should be 
replaced by hi. 

h. Average Interior Eoof Overpressure P. . Here again it is 
necessary to place the structure being considered into one of the tvo cate- 
gories established in the previous paragraphs due to the difference in be- 
havior of the inside shock. For the smaller interiors the average over- 
pressure on the inside of the roof vill build up in time L./U. to the 
initial peak overpressure of the inside shock wave, at which time reflec- 
tion of the interior shock from the rear face vill take place. This re- 
flected shock then travels toward the front vail spreading higher pres- 
sures along the underside of the roof. The reflected shock front is 
destroyed by the vortices at the front face as well as by flow out through 
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openings and no fur- 
ther reflections oc- 
cur. The average roof 
overpressure then de- 
"creases from the rear 
face reflected value 
to the equilibrium 

value of P . the in- 
s 

cident air blast 
overpre s sure . Figure 
3.46 displays the 
average interior roof 
overpressure for 
those structures in 
vhich (I^/U^) 
< O.lt . 
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P i-refl 





i 

ICL" 
T 

u. 

3 





Figure 3.46. Average interior roof overpressure vs timesmall 
rectangular structure (L./ 1). ) < 0. It 



Figure 3 AT dis- 
plays the average in- 
terior roof overpres- 
sure curve for those 
structures in vhich 
yu lo >0.11 . In 
these structures the 
average overpressure 
does not build up to 



Time, t 



Figure 3.47. Average interior roof overpressure vs timelarge 
rectangular structure (L./'U. )> O.lt 



P i-refl since the 
shock wa-'e reflected 



from the rear face 
decays before it can 
spread its higher pressures over the whole of the underside of the roof. 

i. Average Net Roof Overpressure P . . The average net over- 
r ne'e 

pressure acting down on the roof P is determined by subtracting the 

r ne~t 

average interior roof overpressure from the average exterior roof 
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Average Interior Roof 
Overpressure, T*j- ro0 f 

Average Exterior Roof 
Overpressure, P r0 of 

Average Ntt Roof 
Overpressure, P r - nt t 




3-nj 

overpressure on a common time basis 
according to the equation 

P r-net * P roof ~ P i~roof ^ 3 ' 22 ^ 
Figure 3A8 illustrates graphically 

the variation of the average net 
roof overpressure with time. 

J* Average Exterior Side 



Wall Overpressure P ( 



side 



The exte- 



rior loading on the side can be 
handled as for closed rectangular 
structures in paragraph 3-08f , ex- 
cept that the quantity h f should 
be replaced by hi. It is assumed 
that the side is without openings; 
however, even if windows are present, the combined action of exterior and 
interior pressures nay prevent their breakage. In any case, the effect of 
openings is considered to be negligible. 



Figure 3.48. Average net roof overpressure vs 
time rectangular structure with openings 



k. Average Interior Side Wall Overpressure P, 



i-side* 



The average 



interior side wall overpressure is handled in exactly the same manner as 
the average interior roof overpressure* It is assumed that the side is 
without openings, as for the exterior side wall overpressures. 

1- Average Net Side Wall 
Overpressure P ,. The average 



net overpressure acting inward on 



P tldt -f5,. 



sid* 



the side wall P 



is deter- 



's-net 

mined by subtracting the average in- 
terior side wall overpressure from 
the average exterior side wall over- 
pressure on a common time basis ac- 
cording to the equation 

P s-net * P side ~ P i-side (3-23) 
Figure 3^9 illustrate* graphically 

the variation of the average net 
side wall overpressure with time. 



Average Interior Side 
Wall Overpressure, Pj. tldt 

Average Exterior Side 
Wall Overpressure, P g j dt 



Average Net Side Wall 
Overpressure, 




figure 3.49. Average net side wall over- 
pressure vs time rectangular 
structure with openings 
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|P 3-12 PROCEDURE FOR COMPUTATION OF LOADS OH RECTMTGULAH STRUCTURES WITH 
OPENINGS. Repeat Steps 1, 2, and 3 'under paragraph 3-09. 

a. Average Exterior Front 'Wall Overpressure vs Time. Step 1. De- 
termine h' from the relation 



f *-* A f 

Step 2. Determine h r , the clearing height for the front face, treat- 
ing it as if it were closed. Use h' * h ! , if Step 1 gives h 1 > h 1 . 

Step 3. Determine t r from the relation t r 3hL/c . 

Step k. Use Steps 1, 2, and 6 to 9 of paragraph 3-09a to obtain the 
curve of P front vs time* 

b. Average Interior Front' Wall Overpressure vs Time. Step 1. De- 
termine h' f the T, ... &hted average build-up height for the inside front 
wall. 

Step 2. Determine ? QO > knowing P go and A /A , from figure 
3.kl or figure 3-^-2. 

Step 3. Determine U. , knowing P , from figure 3.9. 

Step k. Determine p iMrefl ^ knowing P t > from figure 3.20. 



Step 5 Determine the curve of P using figure 3.22 or table 3.1. 

s 

Step 6. Plot the curve of P. - , vs time as illustrated in 

i-front 

figures 3-38 or 3-39- 

c. Average Net Front Wall Overpressure vs Time, Step 1. Knowing 



P front and P back f r a se( l uence of "ti^^s -b , t^, t , etc., determine 



" back for each ^ lue of * 

Step 2. Plot the curve of P-, , . vs time as illustrated in 

r-net 

figure 3.hQ. 

d. Average Exterior Back Wall Overpressure vs Time. Step 1. De- 
termine hJ from the relation 

6 h A 



Step 2. Determine h 1 , the clearing height for the rear face, treat- 
ing it as if it were closed. Use hJ - h 1 , if Step 1 gives h* > h ! . 
Step 3. Determine t from the relation t* ^^,1/c . 
Step ^4-. Use Steps 1, 2, and k to 8 of paragraph 3-09b to obtain the 



curve of Fl , vs time. 
back 
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e. Average Interior Back Wall Overpressure vs Time, Step 1. Deter- 
mine h* , the weighted average build-up height for the inside back wall, 

Step 2. Determine P gol j U io > and * > i _ re:fl a s in paragraph 3-121). 

Step 3* Determine t d * L j/ U i o * 

Step k. Determine overpressure at time t - t, L./c from the 
relation 

* P soi + C 1 - A of/ A gf )2 (P i-refl ' P soi> 



Step 5. Determine the curve of P + 0.85q as in Steps 1 and 6 to 8 

s 

of paragraph 3-09a. 

Step 6. Plot the curve of P. , , vs time as illustrated in 
figure 3.^3 or 3* l &' 

f * Average Net Back Wall Overpressure vs Time, Step 1. Knowing 
P i-ba k and P ba k for a se< l uence of times ^-^9 t > t , etc., determine 

for ea <* llie of * 



_ 
Step 2. Plot the curve of E. + vs time as illustrated in 

figure 3-^5* 

g. Average Exterior Roof Overpressure vs Time. Step 1. Use Steps 
1 to 8 of paragraph 3-09e to determine the. curve of P vs time, but 
replacing the quantity h ! appearing in any of these Steps by hi . 

h. Average Interior Roof Overpressure vs Time, Step 1. Determine 

P soi> U io> and P i-refl as in P ara S ra P h 3-12b. 
Step 2. Deteraine t d L /2(J . 

Step 3. Determine the curve of P using figure 3.22 or table 3.1. 

s 

Step 1*.. Plot the curve of P i-roof vs time as illustrated in 
figure 3.46 or 3*V7. 

i Average Net Roof Overpressure vs Time. Step 1. Knowing P 
P i-roof for a se< l uence of times \> tg, t 3 , etc., determine P r . net 

" 5 i-roof for each valMe 5 * 

Step 2. Plot the curve of P, , as time as illustrated in 

u-net 

figure 3.1*8. 

J- Average Exterior Side Wall Overpressure vs Time. Step 1. Use 

the procedure of paragraph 3-09g to determine the curve of P _ vs time, 

side * 
but replacing the quantity h f by h f . 
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k. Average Interior Side Wall Overpressure vs Time. The curve of 

p. 4 , is obtained in the same manner as P 4 __. as presented in para- 
i-side i-roof 7 * * 

graph 3-12h. 

1. Average Net Side Wall Overpressure vs Tiine Step 1. Knowing 

de f r a se ^ uence of times \> *> 2 > t-, etc., determine 

PJ, . -, for each value of t. 
i-side _ 

Step 2. Plot the curve of P , vs time as illustrated in 

"* 

figure 3.^9. 



p . P 4J 
s-net side 



LOADING ON STRUCTURES WITH GABLED ROOFS 



I 



3-13 LOADING OH GABLED ROOFS, a. Average Front Slope Roof Overpressure 
P , Shock Front Parallel to Ridge Line. When the shock wave strikes a 
roof surface, or the earth cover of 
a semiburied structure, pitched as 
in figure 3-50(a) or (b), oblique 
reflection of the shock front occurs, 



The reflected press-ore P 



r-a 



is a 




function of the peak incident over- 
pressure P and the angle of in- 
clination of the roof 0. The ratio 

P /P is plotted in figure 
r-or so 

3.11 as a function of the angle of Figure 3.50. Typical gabled roofs 
incidence a, here equal to 90 - 6. 

For a given angle of inclination of the roof, the ratio of the maxi- 
mum average overpressure to the reflected overpressure is a constant when 
the overpressures are expressed in terms of their deviation from the over- 
pressure in the incident shock vave. This relation is expressed as follows, 
where Kg is independent of the peak incident overpressure 



where 



P 1 
roof 



'sL 



peak average overpressure on front slope of each bay at 

b a 



overpi 

time t t , + LA-U 
' 



length of one gabled roof span 

overpressure in the incident shock wave at time t - t. 
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3-13a 



P = reflected overpressure as obtained from figure 3-H 



0.6 



0.5 



S 



. OL 



0.4 



0.3 



0.2 



O.I 



K~ is plotted as 

C7 



10* 20* 30 40* 50 
Inclination of Roof, 6* 



60 70* 



Figure 3.51. Ratio defining peak average 
roof overpressure vs inclination of roof 



The factor 
a function of 6 in figure 351- 
The peak average overpressure on 
the front roof slope pt oof is 
then computed from the following 
equation 

P + K (P - P ) 
sL 9 r-a sL 

(3.25) 

The average overpressure on 
the front slope of the roof builds 
up linearly from zero at time 
t t , - L/4U to the value given 



P 1 _ 
roof 



by equation (3.25) at time t t d + L AU Q - 
clears to the drag overpressure at time t < 
the relation 



for all times in excess of t 

P 



b roof 

. t d + 



t , + 
a 



where 



average overpressure then 



and is given "by 



(3.26) 



roof 

P 
s 



average overpressure on 
the front roof slope 



o -0.6 
o 



overpressure in the inci- tf -0.4 
dent shock wave 



- time-displacement factor 
L/^U for front slope 
of first bay 

q - dynamic pressure 





-0.2 



o 

o +0.2 
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45* 30* 15* 15* 3O* 

Windward Roof Leeward Roof 



45* 



CL drag coefficient plotted 
in figure 3*52 as a 
function of [22]. 

The local overpressure vs time 
curve at a point on the front slope of the gabled roof surface consists of 



Figure 3.52. Drag coefficient for gabled 
roofs vs slope angle of roof 



an instantaneous overpressure rise to the reflected press-ore P 



r-a 



fol- 



lowed by a linear decrease to the drag overpressure, given by equation 
(3.26), in a clearing time equal to L/ 2U * The time displacement for the 
drag pressure curve is equal to the time required for the shock to travel 
from the front edge of the structure to the point at which it is desired 
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L Is Minimum Length of Roof 
Parallel to Shock Front for 
which Curve Is Applicable. 
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to determine the local overpressure. 

The local overpressure vs time curve at a point on the rear roof 
slope consists simply of the drag overpressure given by equation (3.26)* 
The time displacement t, for the drag pressure curve is equal to the time 
required for the shock to travel from the front edge of the structure to 
the point at which it is desired to determine the local overpressure. 

Figures 3.53 and 3-5^ 
illustrate average overpressure 
vs time curves for the front 
i*oof slope in the first bay and 
In any bay. The first bay time- 
displflcement factor t , is 
given above, and for any other 
bay it is the time required for 
the shock front to travel from 
the front wall of the structure 
to a point halfway across the 
front roof slope of the bay. 




K ^< Q \ 

BlasTT f^ Y^\^ 

"to? JW" jpyr. ._ //v 

L- a ^|L/2 L/2 , 



P'roof 



la. 



Note: td - (a + L/4)/Uo 




Time, t 

Figure 3.53. Average overpressure on windward 

roof slope of first bay of gabled structure vs 

timeridge line parallel to shock front 

The local overpressure 
vs time curve at a point on 
the front slope of the gabled 
roof surface consists of an 
instantaneous overpressure 
rise to the reflected pres- 



o + L/2 o + L 
U U 

Time, t 



Figure 3.54. Average overpressure on windward roof 

slope of any bay of a gabled structure vs time 

ridge line parallel to shock front 



sure P followed by a 
linear decrease to the drag 
overpressure, given by equa- 
tion (3.26), in a clearing 
time equal to L/ 2 U Q - The 
time displacement for the 
drag pressure curve is equal 
to the time required for the 
shock to travel from the 
front edge of the structure 
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to the point at vhich it is desired to determine the local overpressure. 

b. Average Rear Slope Roof Overpressure P roof > Shock Front Parallel 
to Ridge Line. No reflection takes place on the rear roof slope, slope B 
of figure 3.50. The overpressure builds up linearly from zero at time 

^k^ 



t, - L/4U , the time at which 
d ' o* 



NOTE: t d 3L/4U 



Blast 
Wavt 



i i 



L l Mlnlnum Length of 
Roof Parallel to Shock 
Front for which This 
Curve Is Applicable. 




the shock front reaches the roof 
peak, to the overpressure as given 
by equation (3*26) at time t t , 4 

L/lrtJ . For the rear slope of the 
o 

first bay the time-displacement 
factor is t, 3LAtT_. The aver- 



.k-'d - 
2U 



. 
U 



'o+'d 



Time, t 



Figure 3.55. Average overpressure on leeward 

roof slope of first bay of gabled structure vs 

timeridge line parallel to shock front 






Ji / I ^^ 4.0,^WW4. 0.0 V^ ^_ / ~Q- 

b i j , ^^ age overpressure-time curve on this 

/ i i ^ s ^>^ 

surface is illustrated in figure 

3-55- 

The local overpressure vs 
time curve at a point on the rear 
roof slope consists simply of the 
drag overpressure given by equation 

(3.26). The time displacement t, for the drag pressure curve is equal to 
the time required to travel from the front edge of the structure to the 
point at vhich it is desired to determine the local overpressure. 

c. Average Overpressures, Shock Front Perpendicular to Ridge Line. 
The overpressures on all surfaces, for the blast wave traveling along the 
ridge line of the structure, are obtained by the methods outlined in para- 
graph 3-08. For closed structures, the clearing height h f is taken as 
the average height of the front wall, or half the width, whichever is the 
smaller; the build-up height of the back wall is taken as the average 
height of that wall or half the width, whichever is the smaller. The 
clearing and build-up heights for gabled structures with openings are com- 
puted by the procedures in paragraph 3-11. 

3-1^ PROCEDURE FOR COMPUTATION OF LOADS ON GABLED ROOFS. Repeat Steps 1, 
2, and 3 under paragraph 3-09. 

a. Average Front Slope Roof Overpressure in First Bay vs Time, Shock 

Front Parallel to Ridge Line. Step 1. Determine U , knowing P , from 

o so 
figure 3.9. 
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Step 2. Determine t d . 



Step 3. Determine K Q , knowing 0, from figure 3.51. 
Step U. Determine P SL , the overpressure in the incident shock wave 
at time t - t d = lAU Q . Use figure 3-10 to get t Q and figure 3.22 to 

find P T /P 
sL so 

Step 5. Determine ^ r ^ knowing a = 90 - e and P , from 
figure 3.11. 

Step 6. Determine P,l on .pj the average maximum overpressure on the 

roof at time t = t, + L/lrtJ , from the relation 

P f = K (P - P ) + P 
roof 9 r-a sL sL 

Step 7- Determine C-y knowing 6, from figure 3-52. 
Step 8. Determine P + C-A by use of Steps 1, 6, and 7 of para- 
graph 3-09a and the following table. 

time, t t - t a (t - t d )/t o q/q. Q 4 C^ P S /P SQ P g P g + C D C 
(1) (2) (3) 00 (5) (6) (T) (8) (9) 

Column (l) is a sequence of times from t = t , to t = t, + t . 

d do 

Column (!.) is obtained from figure 3-2^ or table 3.2 for the values 
in column (3). 

Column (7) is obtained from figure 3-22 or table 3.1 for the values 
in column (3)* 

Step 9 Plot the curve of P vs time as illustrated in 
figure 3-53- 

b. Average Rear Slope Roof Overpressure in First Bay vs Time, Shock 
Front Parallel to Ridge Line. Step 1. Determine U , knowing P > from 
figure 3.9. 

Step 2. Determine 

Step 3. Determine P + C q. by following Steps 7 and 8 of paragraph 

s u 

O _ 1 /In 

Step h. Plot the curve of P vs time as illustrated in figure 

roof 

3.55- 
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LOADING ON EXPOSED STRUCTURAL MEMBERS 



3-15 LOADING ON EXPOSED STRUCTURAL FRAMEWORKS. The blast wave loading on 
structures with open steel frameworks such as bridge trusses, steel towers, 
or steel industrial buildings following the destruction of frangible walls 
is caused almost entirely by the drag pressure resulting from the high- 
velocity blast wind which follows directly behind the shock front. The 
load due to the unbalanced pressures resulting from reflection and diffrac- 
tion of the shock front can be neglected because of the extremely short 
duration of these effects. A structure of this type consists of structural 
elements whose dimensions are so small that their front face clearing times 
are less than 10 milliseconds. In addition the menfbers are so small that 
they are enveloped by the blast wave in a similar short length of time. 
The total load on the structure due to the drag forces can be obtained by 
computing the load on all the individual elements of the buildings $ front 

and back walls and intermediate ele- 



ORAG COEFFICIENTS, C ,FOR VARIOUS 
STRUCTURAL CROSS -SECTIONS 



Co 
2.0 
1.0 
2.0 
2.0 
2.2 



Direction of Wind 

Shape C Shape 

D 2.0 r~ 

I 2.0 O 

I * 2.0 t 1^ 

HH* 1.8 t -i|r 

L 2.0 t Z3 
~1 1.8 

* Standard or Wide Flange Sections 
t Built up Sections, either Riveted or Welded 



DRAG COEFFICIENTS, (C ) w , FOR WINDWARD (UNSHIELDEt 
TRUSS OF BRIDGES 

(C ) w is given as a function of the solidity ratio, G S A'/A, 
where A is the total area included within the limiting 
boundaries of the truss, A 1 is the actual area of the 
members in a plane normal to the wind direction, 
For single span trusses! 

G<0.25 (C ) w * 1.8 
G>0.25 (C ) w 1.6 
For multiple span or very long span trusses! 



0.00<G<0.20 
0,20<G<0.30 
0.30<G <0.90 
0.90<G < 1.00 



(C ). -2.0 
(C ). 1.8 
(Co). ' -6 
(C ) w 2.0 



Figure 3.56. Drag coefficients for various 
structural shapes and trusses 



ments, and adding them together on a 
common time basis. The total drag 
overpressure on an element can be 
computed from the formula 

Drag overpressure CLq, (327) 
where 

CL * coefficient of drag ob- 
tained from figures 3.56 
and 3.57o This coeffi- 
cient is applied to the 
individual structural 
elements 

q dynamic pressure due to 
the air velocity in the 
incident blast wave at any 
time m pu^/2, and is de- 
termined from figures 323 
and 3.24 

t, L/U , where L is the dis- 
tance from the front por- 
tion of the structure to 
the element for which the 
load is being computed 
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0.4 0.6 

Solidity Ratio, G 

Figure 3.57. Drag coefficients for leeward trusses 
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Time, t 

Figure 3.58. Drag overpressure vs 
time on exposed structural element 



3-16 

Figure 3*58 illustrates the drag 
overpressure -time curve for an exposed 
element located at distance L from 
the front portion of a structure. The 
total load on the element at any time 
is obtained by multiplying the drag 
f * td overpressure (at any time) by the pro- 
jected area of the member transverse to 
the direction of travel of the blast 
wave. The total load on the structure 

at any time is obtained by summing the loads on all the individual members, 

based on time t, referenced to the instant the frontmost elements of the 

structure are struck by the shock wave. 
The overpressure-time curve 

for the front wall supporting ele- 

3 

ments of a structure with frangi- S 
ble wall covering is shown in $ 

o 

figure 3*59* Tbe overpressure on ? 

Q 

the front wall prior to the time 
of failure is determined as pre- 
viously described for the front 
wall of a closed rectangular 



Time, t 



Figure 3.59. Drag overpressure vs time on front 

wall elements supporting frangible wall panels 



structure. The time of failure is actually a very short time for frangible 
elements compared to t # and the load transmitted to the supporting ele- 
ments prior to failure may be neglected* After the front wall fails the 
supporting elements are subjected to the drag overpressures given by equa- 
tion (327). The overpressure-time curve for the back wall supporting 
elements is similar to that for the front wall except all times are dis- 
placed by the time -displacement factor t , L/U , where L is the dis- 
tance between the front and back wall. After failure of the back wall, the 
back wall supporting elements are subjected to the drag overpressures 
given by equation ( 3t 27) . 

3-16 EROCEDURE FOR COMPUTATION OF LOADS ON EXPOSED STRUCTURAL FRAMEWORKS, 
Repeat Steps 1, 2, and 3 of paragraph 3-09. 
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Determine U 9 knoving P } from figure 3.9. 



Step 

Step 5. Determine t^ time for shock wave to travel from reference 
point to member under consideration. 

Step 6. Determine CL from figure 3.56 or figure 3-57. 

Step 7. Determine a, knowing P , from figure 3.23. 

o so 

Step 8. Determine q, knowing q . for a sequence of times from 
figure 3.2k or table 3.2. 

Step 9. Evaluate CLq. for these times. 

Step 10. Plot the curve of net drag overpressure vs time as illus- 
trated in figure 3.58. 

LOADING OF CYLIHDRICAL SURFACES 

3-17 LOADING ON CYLINDRICAL ARCH SURFACES, a. Cylindrical Arch Overpres- 

sure P ., Axis Parallel to Shock Front. 

Figure 3*60a is a definition sketch for 

the notation used in the determination of Wove 



Blast 



the loads on a cylindrical arch surface. 
Structures in this category are illus- 
trated in figure 3*60"b. These consist 




I / X \ 1 

V^^^"//////&srf^^ 




x s -jij- (cos 1 - cos 6) 

Figure 3.60a. Definition sketch for 
cylindrical arch notation 




(b) 




Figure 3.606. Illustration of structures 
with cylindrical arch surfaces 



of exposed arch structures and sur- 
face structures of various shapes 
covered "by an earth fill. If the 
arch or earth-fill surface is not 
truly circular, it is approximated 
by an arc of a circle as illustrated 
by dashed lines in figure 3-60b(c), 
and the loads are computed for the 
equivalent circular surface. Note 
that in some cases a better approxi- 
mation of the ground surfade is ob- 
tained with a gabled roof shape , in 
which case the loads are computed as 



75 



3-17a 



1 July 59 



described in paragraph 3-13 Procedures given in this paragraph determine 
the air "blast loading as a function of time on the cylindrical arch surface 
exposed to the air "blast , whether it "be the surface of the arch on the ex- 
posed structure, or the surface of the rounded fill over the covered struc- 
ture. With this procedure j zero time is the instant at which the incident 
shock front first strikes the cylindrical arch surface at its intersection 
with the horizontal ground surf ace , along the line defined by = B l . 

Figures 3-6l and 3-62 illustrate the variation of the local overpres- 
sure P _ with time on the surface of a cylindrical arch* Each of these 
cyl J 



O 




Time, f 



Figure 3. 61. Local overpressure on front sur- 

face of cylindrical arch vs time axis par<- 

allel to shock front, < 0'<0 < 90 



Time, t 

Figure 3.62. Local overpressure on leeward 

surface of cylindrical arch vs. time-~axi$ 

parallel to shock front, 90 < < 180 

figures is valid only within a certain range of 9 values^ which are de- 
termined by the clearing time peculiar to each* These times and the se- 
quence of overpressures within them are defined as follows: 

For < a < 9 < g^_(figjjrg_3^ll:_ The overpressure at a given 
point rises instantaneously from zero to P at time t = x/U The 
value of P is given in figure 3.11 as a function of a, the angle 
of incidence of the shock wave* For the cylinder of figure 3.60a^ a = 9. 



Clearing of reflection effects occurs at time 



t s t, + t 

d c 



where the 



clearing time t is given "by 



t 

c 



_ 
refl 



(3-28) 



The local overpressure at any point after the clearing of reflection 
effects is given by equation (3-29). P acts normal to the surface of 



the arch. 



cyl 



cyl 



(3-29) 
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where q. is the dynamic pressure incorporating the time-displacement 
factor t^ 

C B is the local drag coefficient^ vhich is a function of Q and is 
obtained from figures 363a* 1>j c, and d and 3.64 as discussed 
"below 

P is the overpressure in the incident "blast -wave at any time 

s t-t d 



o 

u 



O 

O 



a. 

u 

"E 

o 




-1.2 



20 



40 



60 80 

(degrees) 



!00 C 



I20 C 



I40 C 



Figure 3.63a. Local dynamic pressure coefficient for cylindrical 
axis parallel to shock fronts-supersonic Mach numbers (M > 1) 
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1,6 




-S.6 



20 C 



40 ! 



60 80 

8 (degrees) 



!00 C 



120 



140 



Figure 3.63b. Local dynamic pressure coefficient for cylindrical arch axis parallel 
to shock fronthigh subsonic Mach numbers (0.4 < M < 1) 



The distribution of dynaiaic overpressures about an arch is a function' 
of both the Reynolds number E and the Mach number M of the high- 
velocity -wind in the blast wave [U8, 5 1 *-]. ^ e Reynolds number, a dimen- 
sionless quantity of the type ~, is- a parameter -which characterizes the 
relative importance of viscous action in steady nonuniform flow. Since the 
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31 

E 



o 
o 




0.2 



0.4 0.6 

Time Ratio, Il 



0,8 



1,0 



Figure 3.64. Mack number of undisturbed stream vs time ratio for various peak overpressures 

viscosity \JL and the density p appear as a ratio in the Reynolds number, 
it is convenient to treat this ratio of fluid properties as a property 
itself. Thus, the Reynolds number involves only the arch diameter D, a 
velocity u, and the fluid property p/(ju The Mach number, also a 
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dimensionless quantity, is the ratio of the actual velocity to the velocity 
of sound. 

In the following discussion, the distribution of dynamic pressures 
about the circular arch is assumed to be the same as the distribution about 
a circular cylinder having the same diameter. Then the results of experi- 
mental studies conducted on flow past circular cylinders can be utilized 
for this problem. .This procedure will introduce some error but it is be- 
lieved to be insignificant compared to the uncertainty associated with the 
basic data. 

For flow at low Mach numbers, the pressure distribution and the total 
drag force change abruptly for a Reynolds number range between 3(10) &&<! 
5(10) . This transition is very sensitive to the smoothness of incident 
flow, temperature gradients, and roughness of the surfaces. The explana- 
tion advanced for this change in pressure distribution and decrease in drag 
is that the boundary layer changes from laminar to turbulent at this crit- 
ical Reynolds number and causes the separation point to move farther back, 
resulting in the redistribution of dynamic pressure observed. 

As the Mach number pf the flow is increased, this decrease in total 
drag and redistribution of pressure at the critical Reynolds number becomes 
less pronounced until at a Mach number of about OA, it no longer occurs. 
The separation point does not move vhen the boundary layer changes from 
laminar to turbulent and is probably fixed by the presence of shock waves. 
For flows with a free stream Mach number greater than OA, the pressure 
distribution is not influenced by the Reynolds number but it is somewhat 
dependent upon the Mach number, the distribution being different for super- 
sonic flow than for subsonic flow. 

Four pressure distributions are necessary in order that the entire 
ranges of Reynolds and Mach numbers can be considered. The critical values 
of Reynolds and Mach numbers are arbitrarily taken as R * 5(10) , M = Q.kj 
and 1.0. These dynamic pressure coefficients are applicable for steady 
flow conditions, and since the flow about the arch is far from steady, 
their use may introduce some additional error. The order of magnitude of 
the error could be ascertained by comparing the results obtained by this 
procedure with direct measurements in field tests in which unsteady effects 
are present. r> 
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The expression for the local pressure coefficient when R < 5(10)^ is 

Q I j- _ -1 

C D = 172 L' 7 + ' 5 ^ - V(L/D) 2 J 



(3-30) 



C' is the local pressure coefficient for an arch infinitely long 

when R < 5(lO) 5 (figure 3-63d) 
L is the length of the axis of the arch 
D is the diameter of the arch 

The Mach number is plotted in figure 3*6^ and the ratio of Reynolds 
number to cylinder diameter is plotted in figure 3.65 for a range of inci 



t - t 



dent peak overpressures as a function of 



i.O 



'd 



I 




10 



10 



10 10 

Reynolds Number Ratio, R e /D 



Figure 3.65. Reynolds number ratio of undisturbed stream vs time 
ratio for various peak overpressures 

For 90 < Q < 180 (figure 3*62); The overpressure at any point on 
the leeward side of the cylinder rises instantaneously from zero to a 
finite overpressure at time t - t , = 0, given by 

P cyl - C 1 '* * e / l80 > P so (3 ' 31) 

This initial overpressure clears to the value of P ^ as given by 

equation (3.29) at time t = t , 4- t , where 

Q. C 

(3-32) 



t c = (Wc refl )(0/90) 
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Appendix A, Transonic Drag Pressures, appearing at the end of this 
manual has "been extracted from a report prepared by Amman and Whitney, 
Consulting Engineers, under contract DA 49-129-Eng-319 with the Chief of 
Engineers. 

A method is outlined and charts included for evaluating pressures oa 
earth- covered aboveground structures that have sloped embankments and a 
limited flat earth fill over their top. The method is based on aerodynamic 
principles which have a laboratory background. The results obtained are 
for a flat-top shock which must be converted to a shock that varies with 
time when applied to the analysis of structures to resist atomic weapons. 
The application to shock waves which vary with time can be accomplished in 
a similar manner described above for the equivalent circle method. At the 
present time the method outlined in appendix A has not been sufficiently 
developed to be set forth as a preferred method of analysis. Supplemental 
data will be available at a future time. 

b. Average End Wall Overpressure P .,, Axis Parallel to Shock 

Front. The average overpressure 

P , on the ends of a cylinder 
end 

oriented as in figure 360a is given 
in figure 366* It is assumed that 
the variation of incident overpres- 
sure in the time required for the 
shock wave to travel across the cyl- 
inder is linear and that the average 
overpressure on the ends is equal to 
the overpressure in the incident 



NOTE: t<j Dcos0'/2U 




0.85q 



Time, t 

Figure 3.66. Average overpressure on closed 

ends of cylindrical arch vs time-axis 

parallel to shock front 



P with the time-displacement factor t , = D cos f /2U . 



blast wave 

c. Local Cylindrical Arch Overpressure, Axis Perpendicular to Shock 
Front. The local overpressure at any point on the periphery of the cylin- 
drical arch oriented with its axis perpendicular to the shock front is the 

overpressure in the incident blast wave P with a time- displacement fac- 

s 

tor t d = L T /U where L T is the distance from the front end to the point 
at which the overpressure is desired. 
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d. Average End Wall Overpressure, Axis Perpendicular to Shock Front * 
The average overpressure on the front and back ends of a cylindrical arch 
oriented with its axis perpendicular to the shock front can "be computed as 
outlined in paragraphs 3-08a and 3~08b The clearing height for the front 
wall and the "build-up height for the "back wall are equal to h as shown 
in figure 3.60a* 

3-18 PROCEDURE FOR COMPUTATION OF LOADS ON GYLMDRICAL ARCHES. Repeat 
Steps 1, 2, and 3 of paragraph 3-09 . 

a. Local Cylindrical Arch Overpressure vs Time, Axis Parallel to 

Shock Front, Step 1. Determine P s knowing a = 9 and P , from 

" ~ ~ ....... r " QJ so 

figure 3.11 for < 90. Initial P = (1.5 - 6/l8o)P for Q > 90. 

cy.L so 



Step 2. Determine c -- , knowing P , from figure 3*21 

re i -L so 

= ^ h / c refl for e < 9 or * c 



Step 3. Determine t = ^ h / c for e < 9 or * = 



(0/90) for >90. 



Step k. Determine U , knowing P , from figure 3.9* 

o so 

Step 5* Determine t, = x/U * 

t - t 
Step 6. Determine 7 - at which M = 1 from figure 3.64, 

"C 

o 
Step 7 Determine CL for times less than this value from figure 

. t - t, 

Step 8. Determine T - at which M = 04 from figure 3.64. 

o 
Step 9 Determine CL for times less than this value "but greater 

than the time determined in Step 6 from figure 

t - t s 

Step 10. Determine T - at which R = 5(lO) p from figure 3*65. 

"C 

o 

Step 11 * Determine CL for times less than this value but greater 
than the time determined in Step 8 from figure 363c, 

Step 12. Determine CL for times in excess of the value determined 



in Step 10 from equation (330) and figure 3.63d. 

Step 13. 
paragraph 3-l4a. 

Step 14. 
3.61 and 3.62. 



Step 13. Determine P + CLq by following the procedure of Step 8^ 

S 13 

Step Ik. Plot the curve of P , vs time as illustrated in figures 
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b, Avemg^^^^^ ^cis Parallel to Shock 

Front. Step 1. Deteimine t d = D cos f /2 u o 

Step 2. Determine values of P for a sequence of times ranging 

s 

from t - t = D cos ! /2U to t - t, = t from table 3,1 or figure 3.22. 
d o do 

Step 3. Plot the curves of P , vs time as illustrated in figure 
3.66. 

c * Local Cyljgdjlc^Arch_0^rpressure vs Time, Axis Perpendicular 
to Shock Front. Step 1. Determine t^, the time required for the shock 
front to travel from the front end of the arch to the point under con- 
sideration* 

Step 2. Determine P for a sequence of times^ with a time- 

s 

displacement factor t . 

Step 3. Plot the curve of P . = P . 

* J cyl s 

d. Average .JBndJHa^_J)v^i2i^sgure vs Time^ Axis Perpendicular to 
Shock Front. Step 1. Determine the clearing height .h ! = h (see figure 
3-63a). 

Step 2* Plot the curves of P , and P, . vs time using 

* front back 

methods of paragraphs 3~09a and b^ respectively* 

LOADING ON DOMES 



3-19 LOADING ON SPHERICAL DOME^SUWA.CES^ The notation convention used to 
designate a point on the surface of a dome is illustrated in figure 3-67- 
With the notation adopted^ any point on the periphery of a spherical dome 
is located by the two angles and 0* Denoting by p any point on the 
sphere, is the angle between the horizontal diameter parallel to the 
direction of travel of the shock and the radius Op joining the surface 
point to the geometrical center of the sphere. The elevation of the point 
p above a horizontal plane through the center of the sphere is R sin d 
sin $? where $f is the angle between the horizontal plane through the 
center of the sphere and the inclined plane containing points 0, 0% and p. 

The variation with time of the local overpressure P, normal to 

dome 

the surface of the dome is similar to that for a cylinder as shown in 
figures 36l and 362o The overpressure at a given point rises 

86 
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at time t, = XU 
- R cos 6 with d being the dis- 
tance across the dome at its "base and 
U is the shock front velocity* The 
angle of incidence at the impinging 
shock wave is equal to For 
values of d less than 90 the re- 
flected overpressure is obtained from 
figure 3*Hj for greater than 90 
the maximum overpressure is given by 
the relation 



View 3-8 



dome 



(3.33) 



The initial maximum overpres- 
sures are cleared in time t as 

c 

given by equation (33^)^ where h 
is the height of the dome above its 
base. 




Blast 
Wav 



Plan Vi@w 



t = 

c 



'refl 



for 6 < 90 



(3.3*) 



Figure 3.67. Definition sketch for 
spherical dome notation 



t = 
c 



for > 90 

The local overpressure at any point after the clearing of reflected 
overpressures is given by equation (335) 



dome 



4- P 



(3*35) 



where 



q. = unit drag pressure incorporating the time-displacement factor t , 

CL = local drag coefficient plotted as a function of in figures 
3-68 and 3,69 

The ratio of the Reynolds number of the air blast wind to the diam- 
eter of the sphere is plotted in figure 3.65 for various incident overpres- 
sures as a function of (t - t,)/t . At the Reynolds number R = 10 , the 

d /7 o * e ' 

value of the drag coefficient changes abruptly. The ratio of (t - t,)/t 



associated with this value of 



R determines the applicability of either 
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figure 3 068 or 3.69 to obtain C D for use in equation (335). 

The loading in domes calculated "by the methods presented in this 
paragraph can be in error for large overpressures* Use of the present 
method should not be made for pressures greater than 15 psi. 
3-20 PROCEDURE FOR COMPUTATION OF LOADS ON SPHERICAL DOMES, Repeat Steps 
1, 2, and 3 of paragraph 3-09* 

Step k. Determine P , knowing a = and P . from figure 3-H 

rQ/ so 

for G < 90. Detemine maximum P. = (1.5 - 0/l8o)P o ^ for Q > 90 . 

dome so 

Step 5 Detemine c -, , knowing P , from figure 3*21. 

rexj- so 

Step 6. Determine t = 3 h / c refl for ^ < 90 or t Q = (3h/c ) 
(0/90) for e > 90. 

Step 7. Determine U , knowing P , from figure 3*9, 

o so 

Step 8* Determine t^ = x/ff . 

Step 9. Determine (t - t,)/t at which R /D = ICr/D, knowing 

do e 

P * from figure 3.65. 
so 

Step 10. Determine CL for times less than the value of t - t, 
obtained in Step 9 from figure 3^68; for times greater than this value 
from figure 3,69* 

Step 11. Determine P + CLq by following the procedure of Step 8, 

S D 



paragraph 3- 
Step 1 
as illustrated in figures 3*6l and 3*62* 



Step 12. Plot the curve of P, vs time for the point selected, 

dome 



LOADING ON BURIED STEUCTUBES 

3-21 LOADING ON BURIED STRUCTURES. The method for deteimination of earth 
pressure loading on buried structures is more approximate than that for the 
loading on surface structures. The soil has an appreciable mass which af- 
fects the magnitude of the overpressures acting on a buried structure be- 
cause of the motion of the structural elements which are in contact with 
the soil* This is handled by an empirical procedure in which the load is 
computed independently of the motion of the structural element. However, 
the mass of the element in contact with the soil is increased by the amount 
of the masj of the soil adjacent to the structure, but of a thickness not 
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to exceed the span of the element. Other factors which complicate the 
situation are n arching action" and wide variations in the density, moisture 
content, and seismic velocity of the soil with depth and location which 
affect the transmission of pressure waves through the ground. The treat- 
ment below, which neglects "dynamic arching," is largely substantiated "by 
the conclusions of reference [46]. 

An air blast wave impinging upon the surface of the ground induces 
underground overpressures of a magnitude approximately equal to the blast 
wave overpressure at shallow depths of 5 to 20 ft [30] in the ground below 
which cause earth overpressures on the sides, roof, and bottom of buried 
structures. The pressure wave thus induced in the ground is propagated at 

the velocity C , the seismic wave velocity in ground (see table 3-3) If 
s 

the ground surface above the buried structures is located in the region of 

Table 3.3. Soil Factors 



Soil Type 


Soil 
Weight 

(lb/ft 3 ) 


C , Seismic 
s 

Velocity 
(fps) 


k, Soil Pres- 
sure Factor 
(Eg. (3-40)) 






Min 


Max 


Min 


Max 


!Bop soil (light dry) 


85 


600 


900 


262 


590 


Top soil (moist, loamy silt) 


95 


1,000 


1,300 


812 


1,370 


Top soil (clayey) 


100 


1,300 


2,000 


1,420 


3,370 


Top soil (semiconsolidated 












sandy clay) 


110 


1,250 


2,150 


1,510 


4,150 


Wet loam 


105 


_ 


2,500 


. 


5,600 


Clay (dense wet, depending 












on depth) 


115 


3,000 


5,900 


8,850 


34,100 


Rubble or gravel 


. 


1,970 


2,600 


6,4oo 


n, 100 


Cemented sand 


14 5 


2,800 


3,200 


9,700 


12,600 


Water -saturated sand 


125 


- 


4,600 


- 


22, 500 


Sand 


3*5 


4,600 


8,400 


26,200 


87,000 


Sand clay 


110 


3,200 


3,800 


10, 000 


13,900 


Cemented sand clay 


145 


3,800 


4,200 


17,800 


21,700 


Clay, clayey sandstone 


150 


- 


5,900 


- 


45,000 


Loose rock talus 


130 


1,250 


2,500 


1,750 


7,000 


Weather -fractured rock 


160 


1,500 


10,000 


3,100 


140,000 


Weather -fractured shale 


150 


7,000 


11,000 


63,000 


156,000 


Weather -fractured sandstone 


160 


4,250 


9,000 


23,500 


116,000 


Granite (slightly seamed) 


170 


- 


10, 500 




io,ooo 


Limestone (massive) 


170 


16,400 


20,200 


390,000 


590,000 
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regular reflection, the air "blast wave will intersect the ground at an 
acute angle, which is assumed equal to zero for computational purposes. 
The induced ground pressure wave front is then horizontal, and moves down- 
ward through the soil and past the structure with a vertical velocity C . 

s 

On the other hand, if the ground surface above the structure is lo- 
cated in the region of Mach reflection, the air shock wave travels parallel 
to the surface of the ground, and the induced ground pressure wave is 
assumed to behave as illustrated in 



figure 3-70. The leading portion of 
the ground pressure wave is inclined 
downward at the angle 7, whose sine 

is C /U , with the surface of the 
s' o' 

ground from the intersection of the 
air shock front with the ground sur- 
face. Although in many cases C 

s 

will exceed U , it is recommended 
for computation that the angle of 



Air Blast Wave 

P. 

M I 1 I i 1 i U 



^Ground Surface 



r sin-' c t /u 

(Assumed Not 
to Exceed 90) 




KP, - 



c s 



Pressure on Underside 
Same as on Roof 



Figure 3.70. Ground pressure wave induced by 
air blast wave in region of Mach reflection 



inclination of the ground pressure wave with the earth's surface not exceed 

90 which corresponds to the case where C = U * As the ground pressure 

s o 

wave impinges on a buried structure, it first loads the upper front corner 
and proceeds horizontally across the roof with a velocity U in the di- 
rection of travel of the air shock wave and simultaneously proceeds verti- 
cally down the front side of the structure with a velocity U X tan j. 
After the ground pressure wave has traveled the length of the roof, it 
proceeds down the rear side again at velocity U x tan y 

When the buried structure is located in either the region of regular 
or Mach reflection, the local overpressure on the roof is approximately 

equal to the overpressure of the air blast wave on the earth's surface P 

s 

since apparently no appreciable reflection occurs on the surfaces of a 
buried structure. The lateral overpressure on the vertical sides of a 
buried structure may be considerably less than the pressures applied at 
the top surface of the soil depending upon the type of soil and the height 

of the water table. A value of 0.15P has been measured in tests at Nevada 

s 

where the soil was a dry siltv clay. Static tests on soils have 
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demonstrated that the lateral at- rest earth pressure varies from 0*4 to 0.5 
of the vertical earth pressure in sandy or granular soils and may become 
as large as the vertical earth pressure for soft clay soils because of a 
phenomena called plastic flow [43, 44, and 45 ] For the type of dynamic 
loading contemplated here, dry clay soils will probably behave similar to 
granular soils because the load will not persist for a sufficient length of 
time to permit the occurrence of plastic flow. However, it is believed 
that wet clay soils would behave like a viscous fluid and transmit lateral 
pressures instantaneously, For design purposes, it is recommended that 
the lateral overpressure on vertical surfaces of buried structures be 
taken as KP , where K can be evaluated as follows: 

K = 025 for dry cohesionless soils 

K = 050 for medium cohesive soils 

K = 0*75 for soft cohesive soils 

K = 1.00 below water table 

a* Local Roof Overpressure P f . Regular Reflection Region: The 
local roof overpressure or overpressure at a given point on the roof is 

equal to the air blast wave overpressure on the ground surface P with 

s 

t, = 0, and where time t is zero when the ground pressure wave strikes 
the roof of the structure . 

Mach Reflection Region: The local roof overpressure is equal to the 
air blast wave overpressure on the ground surface P o with t^ = L ! /U 
where time 

upper front corner of the structure. 
L* is the distance from the point on 
the roof where the pressure is wanted 
to the upper front corner of the 
structure, and U is the air shock 
front velocity* The overpressure- 
time curve is illustrated in 
figuxe 3*71* 

b. Average Roof Overpressure 



s "d ' "o 

t is zero at the instant the ground pressure wave strikes the 



NOTE: 

t d *0 (Regular Refl.) 
t d LVU (Mach Refl.) 



"roof* 



Regular Reflection Region: 



Average roof overpressure is 



Time, t 

Figure 3.71. Local roof overpressure vs time 

for buried structure located in regions of 

regular and Mach reflection 
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10. 

r 



t d L/2U 





Time, t 

Figure 3.72. Average roof overpressure vs time 

for buried structure located in region of 

Mach reflection 

air shock wave. 

c. Local Front and Back Wall 



identical with local roof overpres- 
sure* 

Mach Reflection Region: The 
average roof overpressure is illus- 
trated in figure 3-72. -The over- 
pressure rises linearly from zero 

at time t = to P , the local 

s 

overpressure existing at the center 
of the roof, at time t = L/U , 
where t = L/2U and L = length 

CL O 

of roof in direction of travel of 



Overpressures P. 



front 



and P. 



hack* 



Regular Reflection Region: The local 
front and "back wall overpressures are 

equal to KP . The time-displacement 
s 

factor is t, = h T /U tan 7 where 

d ' o 

h ! is the vertical distance from the 




Time, t 



roof to the point under consideration. Figure 3.73. Local front wall, back wall, and 



side wall overpressure vs time for buried 

structure located in regions of regular 

and Mach reflection 



The overpressure- time curve is illus- 
trated in figure 3-73. 

Mach Reflection Region: The local front and rear wall overpressures 
are equal to KP * The time-displacement factor is t, = h ! /U tan 7 for 

the front face and t, = L/U + h ! /U tan 7 for the rear face where h ! 

d ' o ' o ' 

and L are as defined a"bove. The 
overpressure -time curve is illustrated 
in figure 3*73- 

d. Average Front Wall Qyerpres- 
sure P f . Regular Reflection Re- 
gion: The average front wall over- 
pressure is illustrated in figure 3-7^ 
The overpressure rises linearly from 




Time.t 

Figure 3.74. Average front wall overpressure 

vs time for buried structure located in regions 

of regular and Mach reflection 



zero at time 
94 



t = to KP at 
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t = h/U tan 7, vhere t . = h/2U tan 7 and h is the height of the front 
vail of the structure* 

Mach Reflection Region: The average front vail overpressure is illus- 
trated in figure 37^ The overpressure rises linearly from zero at time 



t = to KP at 



t = h/U tan 
7 o 



where t, 



= h/2U tan 7. 

' o ' 



e. Average Back Wall Overpres- 
sure P, . Regular Reflection Re- 
gion: The average "back vail over- 
pressure is identical vith the aver- 
age front face overpressure* 

Mach Reflection Region: The 
average back vail overpressure is 
illustrated in figure 3*75. The 
overpressure rises linearly from 

zero at time t = L/U to KP 

1 o s 

at time t = L/U + h/U tan 7. 



la. 




Time , t 

Figure 3.75. Average back wall overpressure 
vs time for buried structure located in region 
of Mach reflection 



f . Local Side Wall Overpressure P 



side 



Regular Reflection Region: 



The local side vail overpressure is identical vith the local front and rear 
vail overpressures* 

Mach Reflection Region: The local side vail overpressure is obtained 
in a similar manner as that for the front and rear vail local overpressures. 
The time-displacement factor is t., = L f /U + h ! /U tan 7 vhere L ! is 
the horizontal distance from the upper front corner of the side vail and 
h 1 is the vertical distance from the upper -front corner of the side vail 
to the point at vhich it is desired to obtain the overpressure. The local 
side vail overpressure vs time curve is illustrated in figure 373 

g. Average Side Wall Overpressure P , , . Regular Reflection Re- 
gion: The average side vail overpressure is identical vith the average 
front face overpressure. 

Mach Reflection Region: The average side vail overpressure is illus- 
trated in figure 3*76. The overpressure rises linearly from zero at time 



= ( L / u o + h / u o 



t = to KP g at time t = L/U Q + h/U Q tan 7, vhere 
tan 7)/2. 

h. Overpressure on Underside of Buried Structure. The conclusions 
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reached in reference [52] indicate 
that the overpressures exerted on the 
underside of buried structures with 
an integral floor are of the same 
order of magnitude as the overpres- 
sures on the roof. Therefore, it is 
recommended that in both the regular 



Time, t and Mach reflection regions the "bot- 

Fi gU re3.76. Average side wall overpressure tom surface of b urie d structures with 
vs time for buried structure in region of 

Mach reflection an integral floor be designed to 

support the identical overpressures applied to the roof of the structure as 
given in paragraphs 3~21a and b. 

3-22 PROCEDURE FOR COMPUTATION OF LOADS ON BURIED STRUCTURES. Repeat 
Steps 1, 2, and 3 of paragraph 3-09* 



Step IK Deteimine U , knowing 

o 



P , from figure 3.9. 
so 



Step 5- Determine t, . 

Step 6. Determine P from table 3.1 or figure 3.22. 

s 

Step 7 Plot the curves of overpressure vs time as illustrated in 
figures 3.71 to 3*76. 

RADIATION 

3-23 NUCLEAR RADIATION PHENOMENA. A characteristic of an atomic detona- 
tion is the emission of nuclear radiation. This represents approximately 
15 per cent of the total energy of a typical air burst. It is expressed 
in terms of initial radiation and residual radiation. 

Initial radiation is defined as the nuclear radiation which is de- 
rived directly from the initial fission and fusion reactions of the detona- 
tion. It is delivered within approximately the first half minute following 
detonation and its significant effects are confined to a radius of a few 
miles from the point of detonation. 

Residual radiation applies to the radiation emitted after the first 
minute by the fission products, unfissioned residues, and to a limited ex- 
tent by the bomb case fragments and materials such as dust in the air in 
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which radioactivity has "been induced* Residual radioactivity, in the ag- 
gregate, is of an enduring nature j however, its intensity decreases by a 
process of natural radiation decay* The residual radiation produced by an 
air burst is usually dissipated in the upper atmosphere and does not there- 
fore constitute a hazard on the ground. However, when a nuclear weapon is 
detonated at or near the ground level large amounts of the surface mate- 
rials, such as earth or water, are drawn up into the cloud. Eadioactive 
isotopes formed as described above become associated in various ways with 
these materials which, being considerably heavier than those resulting from 
an air burst, will fall back to the earth's surface in the local areaj the 
heavier particles falling nearer to the point of detonation and the lighter 
materials settling out progressively farther downwind. The material so de- 
posited is known as "fallout" and the areas on which it falls are said to 
be radiologically contaminated* ith megaton yield weapons, fallout may 
extend several hundred miles downwind and cover five to six thousand square 
miles or more. 

Fallout patterns are usually depicted as elongated and cigar shaped. 
Such a pattern would only be found under idealized conditions which seldom, 
if ever, would occur. The existence of complicated wind structures, as 
well as variations of these structures in time and space, may cause extreme 
distortions in the shape of the fallout pattern. In addition, intensities 
within the pattern may be extremely irregular. 

Nuclear radiations consist of neutrons, gamma rays, beta particles, 
and alpha particles. The neutrons, which are subatomic particles of 
neutral charge, are released in the fission and fusion reactions and are 
produced essentially in the first half second after detonation. Gamma 
rays, which are high-energy electromagnetic radiations (like X rays), are 
emitted from the fireball as initial radiation and from fission products 
and from the capture of neutrons in bomb fragments and other materials as 
residual radiation. Beta particles, composed of high-speed electrons 
emitted by the radioactive fission products and alpha particles, identified 
as the nuclei of helium atoms, originate in the unfissioned residues of 
plutonium or uranium. Both the neutrons and gamma rays have a long range 
in air and are very penetrating. Substantial thicknesses of materials of 
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high density are required to reduce their intensity to harmless levels. 
In general, initial gamma radiation is more important than neutron radia- 
tion for large weapons. However, for weapons of small physical size, 
neutron and gamma radiation may be of equal significance, depending upon 
the characteristics of the shielding provided. Beta particles have a very 
short range in air and relatively little penetration capability. Beta 
emitters can damage human tissue when taken into the body, and in addition 
can cause serious burns if they come in contact with the skin. Alpha 
particles have an even shorter range in air and are stopped by almost any 
type of barrier. They are a hazard only if their emitters are taken into 
the body by inhalation, ingestion or, under limited circumstances, through 
breaks in the skin. 

Gamma radiation dosage is measured in terms of a unit called the 
roentgeii r which is a standard measure of the ionization caused by gamma 
rays in their passage through matter and hence of the injury which is 
caused to a body of living organisms. Initial gamma, radiation is a func- 
tion of weapon yield, air density, and the distance of the burst from the 
structure. Figure 3-77 gives the intensity of initial gamma radiation at 
any distance from ground zero for a surface burst. 

Neutron radiation dosage, for the purpose of assessing radiation 
hazard, is measured in rem (roentgen equivalent - mammal) which is the 
amount of energy absorbed in mammalian tissue which is biologically 
equivalent in mammals to one roentgen of gamma or X rays. Neutron radia- 
tion is a function of the weapon yield, weapon design, air density, and the 
distance of the burst from the target. Figure 3-79a gives the intensity of 
neutron radiation at any distance from ground zero for surface burst of 
fission weapons. Figure 3.79b provides similar data for fusion weapons. 

Nuclear radiation can produce a variety of harmful effects in living 
tissue which may be acute or delayed depending upon the total amount of 
radiation absorbed and the period of time within which it is received. 
When radiation is received within a short time (l or 2 days), the effect is 
considered to be acute and, under most conditions, is essentially inde- 
pendent of dose rate. When dosage is received over a long period of time, 
either continuously or in repeated increments, partial recovery takes 
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place. Under these conditions of exposure, larger total doses can be 
tolerated insofar as early effects are concerned* 

Ultimate "body damage resulting from nuclear radiation is a summation 
of the several separate radiation effects involved. Thus, roentgens of 
gamma radiation received must be added to the roentgen equivalent dosage 
for TOR-M (rem) of neutron radiation received in order to obtain the total 
dose. The effects of different amounts of radiation received are shown in 
table 3.4. 

A particular effect is associated with a range of doses rather than a 
specific dose. This does not necessarily mean that any dose within a given 
range will always result in the effect stated but rather that the dose re- 
quired to cause the effect will fall somewhere within this range. 

Table 3.4. Acute Effects of Whole Body Penetrating on Human Beings 



Dose in 

1 Week, 

rem 



Effect 



150 
150-250 

250-350 

350-600 

600 



No acute effects, possible serious long-tern hazard 

Nausea and vomiting within 24 hours, normal incapacitation 
after 2 days 

Nausea and vomiting in under 4 hours. Some mortality will 
occur 2 to 4 weeks. Symptom- free period 48 hours to 2 weeks 

Nausea and vomiting under 2 hours. Mortality certain in 2 to 
4 weeks. Incapacitation prolonged 

Nausea and vomiting almost Immediately. Mortality in 1 week 



In the absence of directives which stipulate otherwise, it is recom- 
mended that protection provided be sufficient to reduce the combined ini- 
tial gamma and neutron dose to not more than 50 rems This will allow for 
possible later additional dosage received from fallout, within structures 
or outside during decontamination operations, or traversal through con- 
taminated areas. 

3-24 NUCLEAR RADIATION SHIELDING. An important consideration in the de- 
sign of a protective structure is the provision of adequate shielding 
against the effects of nuclear radiation. In determining shielding 
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requirements, separate consideration must be given to initial radiation and 
residual radiation. However, ultimate body damage resulting from nuclear 
radiation is a summation of all of the separate radiation effects involved; 
namely, initial gamma, neutron, and residual gamma radiation. 

a. Initial Eadiation Shielding. The effectiveness of initial radia- 
tion shielding is a function of the composition and geometry of the shield- 
ing material, the energy distribution of the radiation at the target, the 
distance from the detonation source, the angle of incidence of the radia- 
tion, and in the case of neutrons, the configuration of the weapon. Ini- 
tial gamma radiation transmission factors for concrete and earth are given 
in figure 3.78. Similar factors for neutrons are given in figure 3.80. 
The curves in these charts are based on the assumption that the radiation 
is essentially point source and is perpendicular to the slab of shielding 
material. The transmission factor is the ratio of the radiation dose re- 
ceived behind the shielding material to the dose which would have been 
received in the absence of the shield. 

The factors which determine the effectiveness of an initial gamma 
shield are density and mass. The transmission factors of materials of 
known density other than concrete and earth may be estimated by interpola- 
tion on a density basis, between the curves shown in figure 378 

Neutron attenuation represents a more complex problem in which 
several different phenomena are involved. First, the fast neutrons must 
be slowed down into the moderately fast range; the moderately fast neutrons 
must then be slowed down into the slow or thermal range; following which 
the thermal neutrons must be absorbed. This absorption process, referred 
to as radiative capture, is accompanied by the emission of relatively high 
energy gamma rays which must be absorbed by the shielding medium. Figure 
3-80 provides a rough approximation of the neutron attenuation capabilities 
of concrete and earth. More reliable data are not presently available but 
studies in this area are in process. Consequently, neutron attenuation 
data should be checked whenever neutron shielding criteria for specific 
projects are required. 

The procedure for calculating the required shielding thickness for 
initial radiation is illustrated by the following example: 
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Figure 3.78. Attenuation of initial gamma radiation 
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Problem: A personnel structure is to be designed that will resist a Iff 

" V" 

50-psi incident overpressure. A 12- in. -thick reinforced concrete roof will 
meet the structural requirement. Determine the earth cover needed to re- 
duce the initial gamma and neutron radiation to an acceptable level of 
50 rem. Consider a 100-KT surface burst at a distance that would produce 
the design overpressure. 

The weapon effects data are: 

100-KT Surface Burst 

50-psi overpressure - distance from ground zero = 0.4 mile 
(figure 3-12*0 

0.4-mile distance - gamma = 70,000 r (figure 3*77) 
0.4-mile distance - neutrons = 100,000 rem (figure 3-79a) 

Procedure: From figures 3-78 and 38Q, the transmission value in 
feet, of concrete and trial thicknesses for the earth cover are found to 
be as follows: 

12 in. of 4 ft of 5 ft of 
Concrete Earth Earth 

Gamma 0.19 0.01 0.003 

Neutrons 0.12 0.003 0.0007 ^ 

Try 4 ft of earth cover. The amount of radiation transmitted is: 

Gamma: 12 in. of concrete, 70,000 X 0.19 = 13,300 r 
Neutrons: 12 in. of concrete, 100,000 rem X 0.12 = 12,000 rem 
Gamma: 4 ft of earth, 13,300 X 0.01 = 133 r 

Neutrons: 4 ft of earth, 12,000 x 0,003 = ^6 rem 

Total transmitted 169 rem 

This is greater than the maximum acceptable level of 50 rem. Try 5 ft of 
earth cover. 

Gamma: 5 ft of earth, 13,300 x 0.003 = 39 r 

Neutrons: 5 ft of earth, 12,000 rem x 0.0007 = 8 rem 

Total transmitted 47 rem 

Since this is less than the maxJuoium allowable of 50 rem, 5 ft of earth 
cover Jls satisfactory. 

"& Fallout Radiation Shielding. The problem of shielding from fall- 
out radiation is different from that of direct radiation. The main differ- 
ences ajre: 
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(1) Eadiation from fallout particles persists for a long period of 
time tut at a delaying rate. 

(2) The radioactive particles are spread over large areas. 

(3) The radiation received at a point on the ground comes partly 
frcm the plane area source and partly from scattered radiations or 
"skyshine." 

(4) The energy of the radiation from fallout is much lover than the 
initial gamma radiation. 

Due to the above characteristics of fallout radiation, the calcula- 
tion of shielding thicknesses is more complex than for direct radiation and 
must take into account the time factor and the height, shape, and size of 
the protected structure together with the location of the point of interest 
relative to the surrounding structure. 

Data and design criteria are contained in references [36^ and Q*-9] A 
detailed method for evaluating the protection afforded by existing "buildings 
against fallout radiation is given in reference [50]* This method can "be 
adapted for use in the design of shielding for protected areas in new 
structures. 

3-25 THERMAL RADIATION, a. Phenomena. The temperatures in the fire ball 
D* an atomic bomb detonation are very high (approximately 330 billion calo- 
ries per kiloton), resulting in a large proportion of the energy being 
emitted as thermal radiation. The thermal energy reaching any particular 
point is a function of the distance from the detonation, the yield of the 
weapons, and the scattering and absorption caused by the atmosphere. The 
mount of atmospheric attenuation is a function of the visibility, or the 
lorizontal distance at which large dark objects can be seen against the sky 
i*b the horizon. The principal characteristics of thermal radiation are 
itaat it: 

(1) Travels with the speed of light (186,000 miles per second). 

(2) Travels in a straight line. 

(3) Has very little penetrating power. 

(4) Can be easily absorbed or attenuated. 

(5) Can be scattered. 

(6) Can be reflected. 
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In view of its speed of travel, the time of arrival at a target is almost 
instantaneous with respect to the detonation* The integrated total thermal. 
energy delivered per unit area is plotted in figures 3.81 and 3.82 as a 
function of distance from the point of burst for a clear atmosphere 
(visibility 10 miles)* Intensity of thermal radiation falling on any sur- 
face is reduced in proportion to the cosine of the angle that the surface 
makes with one which would be perpendicular to the incidence of radiation. 
Thus, the intensity on a surface ^5 degrees from the perpendicular would 
be about 70 per cent of that on a comparable perpendicular surface* 

* Scaling. If two atomic bombs with total energy yields of W, 
and W^ are exploded, the proportion of the total energy emitted as 
thermal radiation is approximately the same for each weapon, therefore, 
the thermal radiation energy is proportional to the yield* At a given dis- 
tance from the detonation the total amount of radiant energy on a -unit area 
is also proportional to the yield and is shown by the expression 



where Q. and Q are the radiant energy on a unit area for the two 
weapons 

c* Effect on Materials* A surface exposed to thermal radiation will 
absorb part and reflect part of this radiation* The amount of the absorp- 
tion depends on the material and especially upon its color, as the darker 
colored materials will absorb a much larger proportion than light colored 
materials* The absorption causes the temperature to rise, and damaging 
effects may result* 

The most important physical effects of the high temperatures due to 
the absorption of thermal radiation are^ of course, charring or ignition 
of combustible materials and the burning of skin* It is very difficult to 
establish definite conditions under which these effects will or will not 
occur, and their extentj however, for convenience, it is found necessary 
to use total energy per unit area criterion* 

Table 35 gives data relative to critical heat energies which produce 
various effects on construction material, surrounding foliage, and exposed 
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Table 3.5. Critical Thermal Energies, Construction Materials 



Material 


Type 


Damage 


Energy, 
cal/sq cm 


100 KT 


10 MT 


Glass 


Ordinary window 


Highly resistant to heat 










"but "brittle and will 










break at low "blast pres- 










sures (see reinforced 










glass) 






Glass re info reed 


Bakelite, cellu- 


Surface melts and darkens 


70 


126 


with plastic 


lose acetate, lu- 








( shatterproof) 


cite, Plexiglas 








Plastics 


USAF window 


Dense smoking and flames 


430 


750 




1/2 in. thick 








Sand 


Siliceous 


Explodes (popcorns) 


19 


35 


Sandbags 


Cotton canvas 


Chars, bags fail 


18 


32 




filled 








Fiberboard 


Corrugated and/or 


Flames during explosion 


10 


18 




laminated 








Plywood 


Douglas fir 


Flames during explosion 


16 


20 


Wood 


Unpainted planks, 


Deep charring 


18 


32 




Douglas fir 








Siding steel, 


28 gage + 


Not affected 
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corrugated 










Siding- galbe st os 


22 gage 


Flames during explosion 


14 


28 


Electric wire, 


12-2-600v 


Flames during explosion 


13 


26 


rub"ber insulated 










BX armored wire 


600v 


Not affected up to 


107 




Romax nonmetallic 


19-2-600v 


Charring 


8 


16 


Roll roofing 


Mineral surface 


Flames during explosion 


ho 


71 


Roll roofing 


Smooth surface 


Flames during explosion 


16 


20 


Grass 


Dry 


Flames during explosion 


7 


16 


Leaves 


Dry fallen 


Flames during explosion 


5 


10 
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3-26 



elements of -structures It will "be noted that more heat energy is required 
from the larger weapon than the smaller one to produce the same effect. 
This is due to the length of the delivery time, which means that in order 
to produce the same thermal effect in a given material, the total amount of 
theimal energy (per unit area) received must be larger for a nuclear explo- 
sion of high yield than for one of lower yield, because the total energy is 
delivered over a longer period of time* 

Although a serious hazard, injury and damage from thennal radiation 
can be more easily avoided than can some of the other damaging effects of 
atomic weapons* Shelter behind any object which will prevent direct ex- 
posure to rays from the detonation is sufficient protection. A wall or any 
other screening object, even clothing, will provide a shield from thermal 
radiation* 



CRATERING 

3-26 CEATERIMG PHENOMENA, When an atomic bomb is detonated near the sur- 
face of the ground^ a crater is formed in the ground under the explosion 
as shown in figure 3-83* Two zones of disturbance in the soil around the 
crater are distinguishable . One, the rupture zone wherein the soil has 
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Figure 3.83. Typical crater 
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been violently ruptured, the other, the plastic zone -wherein the soil has 
been permanently deformed but without visible rupture. If the material 
consists of rock there will be a rupture zone but little or no plastic 
zone. Some of the material thrown out of the crater falls back to form 
a lip around the edge. The height of the lip and the extent of the rupture 
and plastic zones are functions of the crater diameter and depths^ as indi- 
cated in figure 3.83* 

The limits of the rupture zone represent the closest distance to a 
cratering explosion at which it is feasible to construct a resistant struc- 
ture. Due to the detrimental effects of the material thrown out of the 
crater, the high ground accelerations and displacements, and the high in- 
tensity of the blast pressures and radioactivity, it is not considered 
practicable to design protected structures to resist the effects within 
the radius of the crater lip or within approximately twice the crater 
radius. At distances greater than this, the pressure exerted on buried 
structures by the directly transmitted ground shock is less than the earth 
pressure caused by the air blast wave. Although the directly transmitted 
ground, shock is initially of higher magnitude than the air blast wave, it 
s attenuated at a greater rate than the air blast* 

The depth and diameter of the crater depend upon the characteristics 
of the soil and upon the energy yield of the explosion as given by 
figure 38U. 
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APPENDIX A 
TRANSONIC DRAG PRESSURES 

A-01 33KERODUCTION This appendix is concerned with the calculation of 
nuclear blast pressures on aboveground earth-covered structures^ such as 
shown in figure A,l^ which will be applicable to pressure levels above 
25 psig. 

Blast loading 

on structures may 

FLOW^ 

2 



be divided into two 
stages, i.e the 
diffraction phase 
and the enveloped 
phase. When the 
pressure at the 
rear of the surface 
attains the value 
of the pressure in 
the blast wave^ the 
diffraction process 
may be considered 
to have terminated^ 
and subsequently^ 
steady state condi- 




(a) Case 1 

t 

i 

I6 8 



-FLOW 




56' 



i 




^STRUCTURE D 
56' 



(b) Case 2 

Figure A.I. Aboveground earth-covered structures 



tions may be assumed to exist until the pressures have returned to the am- 
bient value prevailing prior to the arrival of the blast wave* The dif- 
fraction phase pressures on the surface of the -earth cover may be computed 
by methods set forth in EM 1110- 345-413 and EM 1110-345-420; however^ little 
information is available concerning the enveloped phase pressures* The en- 
veloped phase pressures consist of the incident pressures plus the dynamic 
pressures. This appendix presents a solution for obtaining the enveloped 
phase pressures on wedge-shaped earth covers similar to those shown in 
figure A.I* The pressures transmitted through the earth to the structure 
may be computed by methods described in EM 1110-345-413 and EM 1110-345-420 
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A-02 



or "by use of a Mohr's Circle Solution. The results are compared -with the 
equivalent circle solution given in these manuals. 

The flows considered are those "behind the shock waves whose charac- 
teristics are given in table A.I. 

Table A.I. Shock Wave Characteristics, 25 psig to 200 psig 



Shock pressure, psig 


25 


50 


100 


200 


Dynamic pressure, psig 


12.2 


41 


123 


330 


Shock velocity, ft /sec 


1756 


2216 


2927 


3985 


Sound velocity "behind shock, ft/sec 


1309 


1447 


1681 


2064 


Particle velocity behind shock, ft/sec 


868 


1375 


2082 


3059 


Temperature at shock, deg Kelvin 

7 / 3 
Density at shock, p x 10 Ib/in. 


379 
871 


467 
1161 


636 
1526 


968 
1895 


Mach number "behind shock 


0.664 


0.951 


1.239 


1.481 


Stagnation pressure "behind shock, psia 


53-5 


116 


294 


766 



The flows and the structures are considered two dimensional. This is 
conservative with respect to the pressures on the structure. The effect of 
the structures being finite in the direction perpendicular to the plane of 
figure A.I is to reduce somewhat the pressures near the ends of the struc- 
ture as -compared to the pressures computed in the following paragraphs. 
A-02 GENERAL FEATURES OF THE FLOW. The general features of the flow de- 
pend on the Mach number of the given conditions behind the original shock 
wave. The Mach number is the ratio of the particle velocity to the veloc- 
ity of sound. When all parts of a flow have a Mach number less than one, 
the flow is called subsonic. If the Mach number is everywhere greater than 
one, the flow is called supersonic. If the Mach number is less than one in 
some places and greater than one in other parts of the same flow, it is 
called transonic. All of the combinations of initial flows and the 
geometries given fall in the transonic range for which both experiments 
and theory are, unfortunately, relatively scarce and complex. 

In the subsonic case the usual concept of pressure coefficients Cp 
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to be multiplied by the dynamic pressure l/2pv are applicable. Coeffi- 
cients given for incompressible flow may be used with less than 1 percent 
error for Mach numbers from zero to O.lo This covers the range intended 
for most wind load coefficients* From Mach numbers 0.1 to 0.5 the incom- 
pressible pressure coefficients may be used if first multiplied by 



This is the Prandtl-Glauert rule* (See reference [A. 2], P !39)* For 
larger Mach numbers, more elaborate methods must be used. 

In the supersonic case the flow pattern is illustrated in figure A, 2. 



SHOCK FRONT 




Figure A.2. Supersonic flow 

On hitting the corner at A, the flow forms a shock wave which turns the on- 
coming flow abruptly, so it is parallel to AB. This causes a pressure 
rise. At B, the flow goes through a Prandtl-Meyer expansion fan which 
again turns the flow, this time with a pressure reduction* At C there is 
.another expansion fan and then a compression shock wave at D returning the 
flow to its original direction. The changes of pressure, density, tempera- 
ture, etc., at each shock and expansion can be conveniently figured using 
figures A. 3 through A. 5* The pressure is relatively constant on each of 
the faces AB, BC, and CD, being of course greatest on AB. The shock wave 
at A is called an attached shock wave because it starts right at A. As the 
Mach number decreases toward one, there will be a certain Mach number 
(figure A. 6) at which the attached shock will no longer be possible for the 
given angle 0. Below this Mach number the shock wave is detached as shown 
in figure A. 7 and the flow will be in the transonic range. 

In the transonic flow shown in figure A. 7, the flow is supersonic 
throughout except for the embedded subsonic region EFBAE. The initial 
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I 



MACH NUMBER (M) 
Figure A.3. Ratio of supersonic stagnation pressure to free stream pressure vs Mach number 
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PRANDTL- MEYER FLOW SOLUTION 
< = >/6 Arc ton V -~^- - Arc tan -/M^l ( Ref. CA.U) 
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Figure A.4. Deflection angle required to expand to a given Mach number 
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Figure A. 5. Stagnation pressure ratio across a normal plane shock (psia) 
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flow is assumed supersonic, 
i.e. M > 1. (Reference 

[A.7]0 

As the initial Mach 
number decreases to 1, the 
detached shock in figure 
A. 7 forms further ahead of 
point A. When the Mach 
number falls below 1, the 
typical situation in the 
transonic range is as 
shown in figure A. 8 with 
one or more supersonic regions embedded in the subsonic flow. The condi- 
tions shown apply to Mach numbers not too close to one. As the Mach num- 
ber increases and approaches one, the supersonic regions grow and coa- 
lesce, finally covering the structure from B to D. In this condition 
there would be a shock wave starting from D as in figure A. 7- 

The flows shown in figures A. 7 and A. 8 have certain features of 
interest in common. In both cases the pressure at A is equal to the stag- 
nation pressure of the subsonic flow upstream of A. At B the local Mach 
number is always unity at the start of the expansion fan. These facts have 



2.0 2.5 3.0 
M-Mach Number 



Figure 4.6. Maximum deflection angle for which the 
shock wave will remain attached (reference \ 
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Figure A.7. Transonic flow (M > 1) 
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A-03a 



SONIC LINE 




V-l EXPANSION 

SHOCK 

( COMPRESSION) 



B ^ C 

SUPERSONIC 



SUBSONIC 



Figure A. 8. Transonic flow (M Q < 1) 

been established experimentally and theoretically (reference [A. 6]). 

Further, it has been shown that the pressure distribution, on the face 
AB is the same for a wide range of transonic flows for a given angle 6 
when the pressure is expressed in terms of the stagnation pressure rather 
than a pressure coefficient (reference [A5]> P 255) "* These facts are the 
basis for the method used below in computing the pressure distribution for 
the cases at hand* 
A-03 METHOD FOR COMPUTATION OF PRESSURES IM THE TRANSONIC RANGE. 

a* Cases In Which The Initial Flow Is Supersonic* The two highest 
pressure shock waves of 100 psi and 200 psi shock pressure yield an initial 
flow which is supersonic as shown in Table A*l* In these cases the condi- 
tions shown in figure A8 apply* 

The first computation is to find the stagnation pressure in the sub- 
sonic region* This may be conveniently done using figures A*3 and A.-5. 
This stagnation pressure at point A, (figure A* 9) is obtained by first 
determining the supersonic downstream pressure f Psi, corresponding to the 
freestream Mach number, M = M^ from figure A*3 (using Y = 1*0, then 
using figure A5 to determine the 
pressure at A, Ps2 For points be- 
tween A and B f for the particular 
angle involved ( = 26 . 6) there are 
fortunately both theory and experi- 
ments available giving the ratio of Figure A. 9. Definition of x and c 
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the pressure at any point x to the stagnation pressure which occurs at A* 
These ratios are plotted in figure 16 of reference A @ 6 and are summarized 
in Table A. 2 

T L? A o, Pressure at x, p f ~ n ^ ^ Q 

Table A.2. 5 !i~ for = 26,6 

Pressure at A, p s 

(From figure 16 of reference \A.6\) 



X 

c 


JL 
ps 


X 

c 


JEL 
ps 





1 


0.6 


0.7^3 


0.1 


0.875 


0.7 


0.723 


0.2 


0.81(0 


0.8 


0.687 


0.3 


0.822 


0.9 


0.655 


OA 


0..793 


1.0 


0.528 


0.5 


0.772 








20 = 53.2 



It will be noted that the experiments and tlieory in most of the refer- 
ences deal with wedges as shown in figure A10 @ The structures considered 

herein are effectively one-half of the 
wedges tested 

The test results and theory avail- 
able are concerned primarily with the 
face AB which is the most heavily loaded 
region and the relatively precise data 
given in Table A* 2 is applicable Un- 
FigureAJO. Typical wedge profile fortunately, equally precise data is 

not available for the faces BC and CD* On the basis of examination of test 
results and of the underlying theory the following procedure is suggested as 
conservative and adequate f though lacking the same precision as the loads on 
AB 

The pressure at B $ on the BC side of the corner may be computed by 
assuming there is a Prandtl -Meyer expansion fan starting from Mach number 
one on the downstream side of the corner as the jiet effect of an over- 
expansion and compensating compression at the corner This pressure may be 
obtained by determining the Mach number f M f on the BC side of the corner 
corresponding to the deflection angle Q = 26@6 from figure A @ k and then 

entering figure A3 with this value to obtain the pressure , P * Along the 

s 

face BC there may be some weak or even quite strong shock waves for which 
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exact data is not available, but their effect will be to increase the 
pressure. A reasonable estimate for the cases being considered is tha^t 
the flow is returned to the pressure of the original supersonic flow up- 
stream of E by the time the corner C is reached. In the absence of exact 
data a linear distribution of pressure between B and C is assumed* 

At C, the pressure on the downstream side may be estimated by 
assuming a Prandtl-Meyer expansion fan starting from the free stream 
Mach number. Along the face BD, it is suggested that the pressure be 
assumed uniform because this region is shielded from the influences that 
can cause pressure rise along BC. The pressure on CD may be obtained 
from figure A.*f and A.5a as follows: (1) enter figure A.*f with free 
stream Mach number to obtain the Prandtl-Meyer angle, 0, (angle tbj^u 
which a supersonic stream is turned to expand from M = 1 to M>*1), (2) 
add (6 to the deflection angle (26.6 for our case) to obtain the ratio 
p/Pt where p is the static pressure on face CD and Pt is the total pressure 
downstream of C. Since there are only isintropic processes between A and 
D, Pt is the same as at A, (3) compute p from the ratio p/Pt. ^ 

b. Cases In Which The Initial Flow is Subsonic The shock waves of ! 
25 and 50 psi in Table A.I have original Mach numbers less than one. In 
these cases, the pressure at A is the stagnation pressure given in Table 
A.I and is obtained from figure A. 3 only by entering the figure with the 
free stream Mach number to obtain the pressure ratio. 

The coefficients given in Table A.2 apply for the region AB, using 
the stagnation pressure as given in Table A. 2 for the angle Q = 26.6 

The pressure on the downstream side of B is computed assuming an 
expansion fan from Mach number one and may be obtained from figure A. 3 
and A.^f for = 26.6 as for the supersonic case. A linear pressure rise 
is assumed from B to C, reaching the free stream pressure at C. 

The downstream pressure at C is estimated according to their being 
a separation of the subsonic flow at the corner and with th$ help of 
experimental data (reference A.9) Simple inviscid theory would indi- 
cate that the flow expands to zero pressure at the corner and thereafter 
resumes free stream velocity. However,, since for these Mach numbers of 
approach to an edge of such 3 large angle as that at C, the flow separates, ^ 
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viscous effects can not be ignored (reference A9) The pertinent data 
in reference A9 (page 28 for station 12, 16) applies to a 30 flap de- 
flection at an angle if attack of the flap of ~9 It would be desira- 
ble to have data for a leading edge flap deflection of 266 and an 



angle of attack of the flap of but the above data is the closest 
available at this time The pressure coefficient f P = (p-p )/q indi- 
cated by the above test data is between -06 and -07t where p = the 
static pressure (psia), P = free stream pressure (psia) and q = the 
dynamic pressure* From another standpoint f if we can apply data from 
a wedge to represent flow on a half wedge following a straight section, 
as we did on the forward face AB, then we can apply data from a sharp 
edged plate at an angle of attack of 26*6 to represent the flow, around 
a sharp corner of 26 6,. with about the same validity* The data on 
pages 18 and 19 (zero flap deflection) of reference A9 show that the 
flow is separated and P is constant over the entire chord for 150, so 
it surely is also for <x = 26*6 The value of P is about 06, and this 
is used for both the 25 and 50 psi cases in this report 
A-04 MJMEBICAL RESULTS FOR THE GIVEN CASES The absolute pressures in 
pounds per square inch (absolute) are shown in figure A. 11 for the given 
shock waves of 25, 50, 100, and 200 psig, respectively* The exact varia- 
tion of the pressure along AB in each case may be found by applying the 
coefficients of table A* 2 to the pressure given for point A (figure A. 9) 
in ach case* For comparison, the pressure distribution obtained using the 
equivalent circle method described in EM 1110-345-413 and EM 1110-345-420 
is shown for the 25-psig level, and the 200-psig level in figures A. 12 and 
A* 13. The values given in figure A. 11 are based upon a "flat top" shock 
but the preceding method is also easily applied to shocks which vary with 
time* 

The pressures given are intended to apply to both of the structures 
shown in figure A.l(a) and figure A.l(b). The difference in the length of 
BC in the two cases will, in fact, make little difference to the pressure 
distribution and the available information does not 'provide any basis for 
computing these variations* It would be desirable to have a solution for 
flow over a double-wedge profile of 53 nose angle, but we have not been 
able to locate any data for this range* 
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PRESSURES, PSI ABSOLUTE. FOR 25PSIG SHOCK WAVE 



64.7 




40 I 



A D 

PRESSURES, PSI ABSOLUTE, FOR 50PSIG SHOCK WAVE 



290 



114.7 




28.0 



A 

PRESSURES, PS! ABSOLUTE, FOR IOOPSIG SHOCK WAVE 
716 




45 3 



A 

PRESSURES, PSI ABSOLUTE, FOR 200PSI6 "SHOCK WAVE 

Figure A.ll. Steady state pressures for the structures in figure A.I 
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EQUIVALENT CIRCLE METHOD (EM 1110-345-413 AND EM 1110-345-420) 
NOTE: ALL PRESSURES ARE PSI ABSOLUTE 



Figure A.12. Comparison of wedge method with equivalent circle method, pressure - 25 psig 
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Figure A. 13. Comparison of wedge method with equivalent circle method^ pressure = 200 psig 
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A-05 REMARKS CONCERNING DOME OR CONICAL- SHAPED STRUCTURES. The above 
analysis has been concerned with structures which are long enough compared 
to the dimension in the direction of flow to be considered two dimensional. 
A dome or conical- shaped structure does not fall in this category and the 
coefficients suggested .above should not be us^d for such cases, although 
this would certainly be conservative as far as the heavily loaded windward 
side is concerned. 

Field measurements on full scale structures of this type were made 
during the 1957 nuclear test series conducted at the Nevada Test Site, 
but were not available at the time of this report. 
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